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ABSTRACT 


When compared to other laser types, the Free Electron Laser (FEE) provides 
optimal beam quality for successful atmospheric propagation. Assuming the 
development and deployment of a mega-watt (MW) class, ground or sea based FEE, this 
thesis investigates several proposed space applications including power beaming to 
satellites, the removal of orbital debris, laser illumination of objects within the solar 
system for scientific study, and interstellar laser illumination for communications. Power 
beaming simulations are conducted within the Satellite Tool Kit (STK) program to 
determine the frequency of accesses between a satellite and one or more ground stations 
for multiple orbital profiles. FEE illumination of orbital debris is modeled to determine 
the thermal effects on a representative aluminum debris particle. FEE illumination of the 
Lunar and Martian surfaces is modeled to determine the relative laser and solar spectral 
intensities at these ranges. FEE illumination at interstellar ranges is modeled and 
discussed to determine our ability to communicate or detect laser communications over 
interstellar ranges. 
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I. INTRODUCTION 


The development of high average power free electron lasers (FELs), with their 
nearly unlimited range and excellent beam quality, offers the potential to revolutionize 
the exploration and utilization of space. This thesis will investigate the utilization of an 
FEE for power beaming to satellites, orbital debris clearance, and illuminating 
extraterrestrial bodies for scientific study or communication. Additionally, recently 
published Chinese forays into anti-satellite warfare suggest other, offensive and 
defensive, military applications of high energy lasers, such as blinding electro-optical 
satellites, engaging ballistic missiles, or damaging enemy satellites in wartime. 
Implementation of such capabilities would violate international law and treaties, and they 
will not be explored here. 

A high average power laser could be used to beam power over large distances to 
Unmanned Aerial Vehicles (UAVs), satellites, or even a future lunar base, but this thesis 
will focus primarily on satellite applications. For satellites, the loss of solar power during 
eclipse often forces the implementation of a restrictive power budget, and additional 
power provided by a high energy laser could affect the operation and design of future 
spacecraft. Utilizing specially designed photovoltaic cells, this method could power 
satellites with much greater efficiency than current solar cell technology. It would 
remain impractical for steady power applications and might only be useful or significant 
for emergent or wartime utilization. 

By combining a high power laser with advanced sensors, a system can be 
deployed that finally matches the range and speed of its sensors allowing for rapid 
illumination of multiple targets. If applied to the growing problem of orbital debris, a 
high-peak power laser could be used to apply a small change in orbital velocity over 
several orbits. By changing the orbital profile to lower the perigee and, therefore, 
increasing atmospheric drag, the laser could greatly decrease the time it takes for debris 
to reenter. Similarly, a high average power FEE could melt and then vaporize some of 
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the debris material, resulting in a smaller and less dangerous partiele and making near- 
Earth spaee safer for satellites and manned missions. 

When applied to extraterrestrial illumination, i.e., the illumination of any non¬ 
terrestrial natural body sueh as the moon, other planets, eomets, asteroids, or even 
interstellar dust, a high average power laser allows for the transmission of long and 
eomplex pulses with good eoherence, leading to improved scientifie observation of these 
bodies over a speeific bandwidth. If applied to the search for extraterrestrial intelligence, 
a high average power, continuous wave laser could be used to send signals to extra-solar 
planets that have the potential to support life. Despite the many difficulties of interstellar 
transmission, the coherence, power spectral density, and directionality of a laser beam 
provides the best opportunity for successfully transmitting a discernable and uniquely 
intelligent signal across interstellar distances. 

Supporting the applications listed above, the first four chapters in this thesis 
review the history and basic theories of laser applications, FELs, laser propagation, and 
orbital mechanics. This background and theory is intended to lay a foundation of 
knowledge and understanding for the analyses performed in subsequent chapters while 
highlighting any specific results or differences that apply to a high power FEE. The thesis 
begins with a basic discussion of laser history in Chapter I, followed by a description of 
the specific qualities, components, and understanding of an FEE in Chapter II. Chapter 
III discusses the propagation of an electromagnetic wave across free space and through a 
turbulent medium, e.g.. Earth’s atmosphere. Chapter IV reviews orbital mechanics, 
describes commonly used orbits, and discusses the perturbative forces that act on an 
orbiting object. 

Chapter V examines the use of a high power FEE for beaming energy to EEO 
satellites. Assuming the successful development and deployment of a high power FEE, 
operational concepts for power beaming are explored and simulated to investigate the 
possible utilization, costs, and benefits of a deployed ground or sea-based infrastructure. 
Chapter VI proposes the use of a high power FEE for the clearance of orbital debris. A 
discussion of the debris types and the motivation for removal is followed by calculations 

of laser energy required to vaporize a significant portion of the debris material. Chapter 
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VII investigates the utilization of a high power PEL for Lunar or Martian illumination 
and possible applications in sending interstellar communications. Building on the 
propagation equations outlined in Chapter III, calculations of the beam spot size, 
bandwidth, intensity, and comparisons to solar illumination at the target are detailed. 
Similar calculations are included for interstellar communications with the goal of 
transmitting a clear and unique signal to extra-solar planets that might harbor life. 

The potential applications of high energy lasers in the not-so-distant future are 
exciting, and each application abounds in interesting scientific theory, engineering 
challenges, and detailed operational considerations. For a naval officer, the militaristic 
applications are especially intriguing due to the laser weapon’s ability to match the speed 
and efficacy of its radar sensor, i.e., the weapon can engage and destroy missiles at nearly 
the same range that the sensor acquires them and the engagement takes approximately 
two seconds instead of ~45 seconds. This thesis performs an initial theoretical and 
operational analysis of a few near-term and several long-term applications of the next 
generations of FELs. 

A, LASER HISTORY 

Despite their deep entrenchment in our popular culture and science fiction, the 
concept of Light Amplification by the Stimulated Emission of Radiation (LASER) and 
directed energy weapons in general are very recent scientific developments. Although 
some reports claim that the Greek scholar Archimedes utilized large mirrors to set fire to 
enemy ships, the first verifiable uses of directed energy were entirely scientific [1]. In 
the late 1800s and early 1900s, the invention of the cathode ray tube and later the 
discovery of the medical utility of X-rays prompted a burst in scientific exploration into 
directed energy. The first true laser wasn’t developed until the mid-20* century, and 
battlefield lasers are only now becoming commonplace. 

Perhaps the most important discovery leading to the development of functional 
maser and laser devices was the principle of stimulated emission proposed by Albert 
Einstein in his 1917 paper “On the Quantum Theory of Radiation.” Prior to Einstein’s 
paper, the scientific community only recognized two fundamental and limited 
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interactions between a photon and an atom in one of multiple available (but discrete) 
energy states. In the first case, “absorption,” the atom was initially in a lower energy 
state when it absorbed a photon with exactly the right amount of energy to move the atom 
into a higher, “excited” state. In the second case, “spontaneous emission,” the atom was 
initially in a higher energy state when it spontaneously dropped to a lower energy state by 
emitting a photon whose energy exactly matched the difference in states. 

Detailing a third alternative in his paper, Einstein postulated that if a photon with 
the appropriate energy arrived when the atom was in the higher energy level, it could 
stimulate the atom to emit another photon of the same energy, and drop to the lower 
energy state. Although pivotal to the eventual development of the laser, this process of 
“stimulated emission” was initially ignored by most scientists because it seemed to be a 
very unusual event. The primary hurdle to the scientific exploration of Einstein’s 
stimulated emission was the concept of “thermodynamic equilibrium.” Eor most of the 
early 1900s, scientists believed that in any population of atoms most physical and 
chemical processes would quickly reach steady-state equilibrium, thereby limiting the 
proportion of atoms that exist in a higher energy state. With most atoms confined to a 
minimum energy (ground) state, the process of absorption would dominate any 
contributions from stimulated emission, and lasing would be impossible. 

Despite these conceptual limitations, scientists of the 1930s learned to utilize 
microwave cavities to generate coherent radiation with an electron beam which prompted 
the rapidly developing radar technologies of World War II. The microwave tube utilized 
a beam of free, non-relativistic electrons and a closed resonant cavity to produce a fairly 
long wavelength electromagnetic beam with good efficiencies. Still, the physical 
limitations of their microwave tubes meant that shorter wavelengths would remain out of 
reach for the foreseeable future, so most researchers focused their efforts on microwave 
wavelengths and technologies. Building on the process of stimulated emission and his 
experience with surplus military microwave equipment, Charles Townes created the first 
“ammonia” maser in 1954 at Columbia University. This maser—short for Microwave 
Amplification by the Stimulated Emission of Radiation—directed a beam of excited 
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ammonia molecules into a cavity resonant at the 24 GHz frequency of the ammonia 
transition to couple the emitted spontaneous radiation with the ammonia molecules and 
thus provide more amplification [2], 

In order to generate coherent radiation at shorter wavelengths, scientists devised a 
way to pump bound electrons in an open optical resonator into excited states within a 
crystal medium, and in 1960, while working at Hughes Research Laboratories, Theodore 
Maiman produced a small flash-lamp pumped ruby laser at 694nm (a wavelength in the 
deep red) [2], Further developments in gas, dye, chemical, and other exotic lasing 
mediums followed until in the 1970s, scientists at Stanford developed the first free 
electron laser (FEL) which utilized a relativistic beam of electrons as a lasing medium, 
thus avoiding the traditional thermal limitations of a static medium [3]. Building on a 
1951 proposal and later demonstrations by Hans Motz in the millimeter wavelengths, 
John M. J. Madey proposed and developed the first FEL, which produced optical waves 
by passing an electron beam and an external carbon-dioxide “seed” laser pulse through a 
series of alternating magnetic fields [2], 

Since their initial development in the 1950s, lasers have become commonplace 
devices in every home, office, and lab; with applications ranging from industrial welding 
to laser hair removal. However, it is often too easy to overlook the scientific and 
technological significance of their invention. In their over half-century of existence, 
lasers have expanded the boundaries of coherent electronics from the millimeter 
wavelength range, using microwave tubes and transistors, to include the submillimeter, 
infrared, visible, ultraviolet, and portions of the X-ray spectrum [4]. Scientists and 
engineers now have access to coherent radiation at frequencies up to a million times 
higher (and wavelengths up to a million times shorter) than previously possible. 

B, BASIC LASER CONCEPTS 

In the broadest sense, a laser could be defined as any device that generates or 
amplifies coherent light through the process of stimulated emission. However, this 
definition requires an equally broad definition of the term “light” due to the wide 
spectrum of wavelengths accessible to modem lasers (from X-ray to microwave). 
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Generally, laser light exhibits remarkable coherence, spectral purity, directionality, and 
intensity, but specifically, there are three essential elements to every laser; a gain 
medium, a pumping process, and an optical feedback process. 

1, Essential Elements 


While ah lasers combine these three basic elements to generate or amplify an 
optical wave, each element or process can be constructed in a variety of ways. The lasing 
medium can be any appropriate collection of atoms, molecules, ions, electrons, or even a 
semi-conductive crystal with multiple quantum energy states. The pumping process must 
have some way of repeatedly exciting the medium into a higher energy state to replenish 
the laser between passes. The optical feedback process must amplify a beam of radiation 
as it passes through the lasing medium over one long pass (as in a laser amplifier) or 
multiple repeated passes (as in a laser oscillator, shown in Figure 1) to maximize the 
interaction between the optical wave and the gain medium. 
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Figure 1. Basic Schematic of a Conventional Laser. From [5] 
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2, Population Inversion 


The stimulated transition proeess within a gain medium ean either absorb or 
amplify an applied optieal signal. In absorption, atoms initially in the lower state are 
stimulated by an applied signal to transition upward to the higher state; in amplification, 
atoms initially in the higher state are stimulated to make the transition downward to the 
lower state. The stimulated transition probability produced by an applied signal is 
exactly the same in both directions. Therefore, although both processes take place 
simultaneously, the net flow of stimulated transitions is always from the energy level 
with a larger population to the energy level with a smaller population. In order to 
produce laser amplification and to overcome the thermodynamic equilibrium described 
above, the pumping process must create a condition called a population inversion, in 
which more atoms are in the higher energy level than are in the lower energy level [4]. 

3. The Pumping Process 

As shown in the upper portion of Figure 1, most conventional lasers utilize the 
difference in quantum energy states within their gain medium to produce an optical beam 
with a narrow wavelength range. By “pumping” electrons into the higher energy states, 
the optical output of both spontaneous and stimulated emission becomes highly 
predictable, but not perfectly so, due to small differences in the allowed energy levels. 
In modern lasers, pumping can be accomplished through many familiar and sometimes 
exotic methods such as chemical reactions, gas discharges, direct electrical pumping, and 
nuclear pumping. In most lasers, however, the pumping process occurs through the use 
of optical pumping methods such as the use of flash lamps or even other lasers. 

4. Coherence 

As an optical signal is applied to the laser’s gain medium, the electromagnetic 
fields within the wave cause a resonant response of the atomic wave functions to the 
applied signal. In a classical sense, each atom of the medium undergoes internal 
oscillations that follow the driving optical signal coherently, matching it precisely in 
amplitude and phase. As these oscillating atoms reradiate, the fields reradiated by the 
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individual atoms combine with the applied fields to produee absorption or amplifieation 
in a manner that is both spatially and speetrally eoherent. 

5, Laser Limitations 

Due to the limitations of the lasing medium, a typieal laser is limited in the range 
of available output wavelengths, the width of the output optieal speetrum, the total output 
power, and the overall system effieieney. The output wavelengths are limited beeause 
there are a limited number of transitions available in a partieular gain medium. The 
width of the output optieal speetrum is affeeted by many slight variations in allowable 
quantum-meehanieal energy states. The optieal output power is limited beeause of the 
thermal stresses within the gain medium. Thermal effeets within the gain medium ean 
deerease the optieal beam quality, and at extremely high powers, the growing thermal 
stresses ean eause damage to the medium or eavity and thereby result in a eatastrophie 
failure of the laser. The overall system effieieney is limited in two ways: first, waste heat 
deposited in the medium must be removed; and seeond, the pumping proeess ean also 
exeite unwanted energy transitions within the medium. These unwanted energy 
transitions produee photons outside of the desired output range, whieh ean eontribute to 
the laser’s thermal loads or degrade the overall beam quality. 
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II. THE FREE ELECTRON LASER 


Although largely differentiated from other elasses of laser in size, funetion, and 
eomplexity, the PEL still eombines the three basie eomponents of a laser to generate and 
amplify eoherent light. In an FEE, the lasing medium is aetually a high energy beam of 
relativistie eleetrons traveling through an undulator at nearly the speed of light. Avoiding 
the thermal eomplieations that plague other laser types in their media, the FEE eapitalizes 
on the prineiple of synehrotron radiation to deliver a high power laser beam limited in 
power primarily by the thermal properties of the laser opties. With this advantage, FEEs 
have the potential to be mueh more powerful than other elasses of laser onee FEE 
teehnologies are fully developed. Eikewise, by freeing themselves of the limited atomic 
transitions within a standard medium, FEEs can be easily tuned over a broad range of 
wavelengths. 

The pumping process in an FEE is provided through the addition of radio 
frequency (RE) accelerating cavities that continually raise the electron beam to 
relativistic energies. These RE cavities and the large cooling systems needed to maintain 
them are the greatest contributors to the added size, weight, and complexity of FEE 
design. The cavities can add significantly to the “wall-plug” efficiency of the complete 
system by recapturing most of the energy of the returning electron beam when the beam 
is recirculated. 

Much like other laser types, the optical feedback process in an FEE can occur in 
both the amplifier and oscillator designs. In the amplifier, a long undulator maximizes 
the interaction time between the electron beam and light generated either by a seed laser 
or by spontaneous emission. In the oscillator, the optical feedback is provided by a 
resonator formed by two mirrors bracketing the undulator in order to maximize the light 
and electron beam interaction. 
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A. 


FEL COMPONENTS 


In terms of the basic laser theory described above, the FEL is categorized as a laser 
because it consists of an electron beam used as a gain medium, a pumping method of 
accelerating these electrons to relativistic energies, and a method to provide optical feedback 
and amplification through prolonged interaction between the electron beam and the optical 
wave. These statements are simplified to connect basic lasing principles to the FEL, but in 
terms of actual components, there are two general FEL designs: the oscillator and the amplifier. 
As shown in Figure 2, the basic components of an oscillator FEL design are: an electron 
source, an electron accelerator, undulator magnets, resonator mirrors, and an electron dump. In 
contrast, the amplifier design has an electron source, an electron accelerator, a longer undulator, 
a seed laser, and an electron dump. As the names imply, the primary difference between 
oscillator and amplifier designs is the removal of the optical resonator; therefore in an 
amplifier, all of the laser’s gain must occur in a single pass through the undulator, and the initial 
optical pulse is often provided by the addition of a seed laser. 



Figure 2. JLAB Energy Recovery LINAC and Oscillator FEL. From [6] 
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1, Injector 


The first and often most limiting eomponent in any PEL is the injeetor, whieh 
generates the eleetrons that later beeome the FEE gain medium. An ideal injeetor would emit 
a narrow and very dense beam of eleetrons for later aeeeleration by the aceelerating cavities. 
Due to the repelling effects of the Coulomb force, larger electron densities tend to increase 
the initial beam emittance, the product of the radius and angular spread of the electron beam. 
FEEs can utilize several types of cathode inside their injector—the photocathode, the 
thermionic cathode, and the field emission cathode are examples. As shown in Figure 3, a 
typical photocathode utilizes the photoelectric effect to liberate electrons from the metallic 
cathode surface by illuminating it with a pulsed laser. Similarly, the thermionic cathode 
imparts kinetic energy to the electrons through a heating element to dislodge electrons from 
the cathode. A field emission cathode relies on the application of an exceptionally high 
external, electric field to produce tunneling and thereby emit electrons from the cathode 
surface. In all of these cases, once the electrons are released from the cathode surface, a high 
voltage field inside the injector’s housing quickly accelerates them away from the cathode 
and toward the accelerating cavities. 
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Figure 3. Typical Photocathode. From [5] 
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2, Accelerating Cavities 


After departing the eleetron gun, the eleetrons are aeeelerated to the required 
relativistie veloeities through the use of radio frequeney (RF) aeeelerating eavities. This 
aeeeleration is necessary to produce an electron beam with enough energy to create laser 
output at the desired wavelength. FELs have operated with linear accelerators (LINACs) 
consisting of either “normal” copper cavities or superconducting cavities, and either 
design has various benefits depending on the FEL’s intended wavelength and output 
power. In both cases, the addition of RE cavities and their associated cooling systems 
adds considerable weight and volume to the overall EEE system design. 

Despite the added cooling requirements, the superconducting accelerator cavity is 
generally preferred for high-power applications [7]. By reducing the resistive losses in 
the cavity wall, the RE energy that would have been lost as heat in a room temperature 
linear accelerator can now be applied to accelerate electrons, thereby improving overall 
cavity efficiency. Eikewise, the higher average-power gradients created by 
superconducting cavities allow for shorter overall accelerating structures. 

The accelerating cavities play an additional role in energy recovery linac (ERL) 
designs by allowing for the recovery of unused kinetic energy in the electron beam. Eor 
EEL designs at optical wavelengths, very high electron beam energy is required, but only 
a small fraction of this energy is actually utilized in the production of light by each 
electron bunch. Therefore, without some mechanism for energy recovery, a large amount 
of initial energy would go to waste in the beam dump. By carefully timing the arrival of 
returning electrons in the accelerator to be approximately 180 degrees out of phase with 
the oscillating RE fields, the electrons will actually decelerate and transfer most of their 
kinetic energy back to the accelerator. This allows for the recovery and reuse of a major 
portion of the original input power and reduces the thermal and radioactive contributions 
at the beam dump. 

The models and picture in Eigure 4 display the construction of a typical 
superconducting RE cavity. The upper figure is a single accelerating cavity with the 
electric fields, magnetic fields, and electron beam graphically depicted. In the center 
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figure, multiple cavities are combined, and the alternating electric fields and electron 
beam are displayed. The bottom figure is a photo of a typical accelerating cavity. 



Figure 4. A Cylindrical RF Cavity with Longitudinal Electric Field for Particle 

Acceleration. From [5] 

3, Cooling Systems 

While a superconducting material exhibits zero loss for direct current (DC) 
applications, alternating current (AC) applications in the RF range produce a loss that 
scales in a roughly exponential manner with temperature [8]. Therefore, in order to 
reduce the resistive losses in the cavity walls, superconducting RF accelerator cavities 
must be cooled from room temperature (-300 K) to nearly 2 K. With such a large 
temperature gradient, these systems commonly utilize multiple insulation layers and 
shields at progressively lower temperatures to obtain the desired cavity temperature. 
Using liquid helium as a refrigerant, the cooling systems move heat away from the 
exterior walls of the accelerating cavities using a standard refrigeration cycle, but the 
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ultimate size and total input power of the cooling system is driven by the size of the 
accelerator cavities and their required cooling rate. 

At the center of the cooling system, a volume of liquid helium surrounds the 
accelerator cavities. As shown in Figure 5, there are two liquid phases of helium—He I, 
which exhibits the properties of a conventional fluid, and He II, which is a “superfluid.” 
Below about 2 K, liquid helium exhibits unusual dynamic properties (thus the description 
as a superfluid) and an extremely large heat capacity. This large heat capacity is what 
makes He II an ideal refrigerant for high power and high frequency applications. 
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Figure 5. Helium Phase Diagram. From [9] 
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At atmospheric pressure, the boiling point of liquid helium is 4.2K, and for lower 
frequency applications, the cavity power density and losses are low enough to use 
“normal” He I for a refrigerant. Working at the natural boiling point allows for a 
simplified cryogenics system, but at frequencies above 800 MHz, the losses and power 
densities generally require more cooling than a He I system can provide. In the case of 
RF applications, the accelerating cavities are immersed in a bath of superfluid He II 
contained by a secondary shield. The secondary shield is generally maintained at 
temperatures in the 30-50 K range by encasing them in a bath of He I. By reducing the 
drastic thermal gradient from the hot accelerator cavities to the room temperature 
housing, the cooling systems gain efficiency and reduce the total input power 
requirements [10]. 

4. Beam Control Systems 

From the moment the electrons are ejected from the cathode until their eventual 
collision with the beam dump, the electron beam radius and angular spread must be closely 
maintained. Inside the undulator, the overlap of each electron bunch with an associated 
optical pulse is essential to the production of light and to the operation of the FEL. 
Capitalizing on the magnetic component of the Lorentz force, beam steering and focusing is 
often achieved with the addition of large quadrupole electromagnets. Other bending and 
focusing systems can use electric or magnetic dipoles to alter the beam’s trajectory. 

As shown in Figure 6, quadrupole magnets consist of four magnets laid out in a 
“cross” pattern. Therefore, in the multipole expansion of the magnetic field, the dipole 
terms cancel and the lowest significant terms in the field equations are quadrupole terms. 
This effect produces a magnetic field that grows rapidly with radial distance from the 
beam axis, thereby providing an innate mechanism for beam focusing. It is impossible to 
simultaneously focus in both the horizontal and vertical planes in a quadrupole, so that a 
set of quadrupole magnets will focus in one plane while defocusing in the orthogonal 
plane in sequence along the beam line. By properly arranging alternating sets of these 
magnets (horizontally focusing then vertically focusing) with appropriate spacing in- 
between them, it is possible to achieve an overall focusing effect on the electron beam. 
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Figure 6. Quadrupole Electromagnetic with Field Lines. From [11] 

5, Undulator 

As shown in Figure 7, the undulator is a collection of strong permanent magnets 
that establish a periodic magnetic field to elicit oscillations in the electron beam. As 
electrons travel along the longitudinal axis, these alternating magnetic fields induce 
oscillations of the electrons in a transverse direction with respect to their original motion. 
Following the process of stimulated emission, the periodic transverse acceleration of the 
electrons in the presence of an optical wave produces coherent amplification of the 
optical wave. In order to ensure that amplification continues, both the amplifier and 
oscillator FEE designs must carefully coordinate the overlap of optical pulses and 
electron bunches. Demonstrating one of the primary advantages of the EEL, it is 
relatively easy to “tune” the output wavelength of an FEE by adjusting the spacing 
between the undulator magnets, commonly referred to as the undulator period. 
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Figure 7. Periodic Longitudinal Undulator. After [12] 


6, Resonator 

Utilized only in the FEL oscillator design, the resonator is an optical cavity created by 
sandwiching the undulator between two reflective mirrors as shown in Figure 8. The distance 
between these mirrors is finely adjusted in order to ensure amplification of the optical wave and 
to prolong the interaction between the optical wave and the electrons. By ensuring the proper 
synchronization of electrons and optical pulses, the optical cavity allows for the amplification 
of the original light pulse. By making one of the mirrors partially transmissive, a fraction of the 
laser light is released on each pass to the output optics and beam director. 
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Figure 8. 


Generic Oscillator FEL Design. Erom [13] 








7. 


Electron Beam Recirculation / Energy Recovery 


While passing through the undulator, the electron beam only contributes a small 
portion of its kinetic energy to the optical wave. There are two primary designs to 
recover this residual kinetic energy and improve the total system efficiency—the Storage 
Ring and the Energy Recovery LIN AC. In an FEE utilizing a storage ring design, 
steering and focusing magnets bend the electron beam back to the accelerating RE 
cavities for re-acceleration back to the required operational energy. Thus, in this case, 
the same electrons will make multiple passes through the undulator. In contrast, an 
energy recovery LINAC ensures that the returning electron beam is properly out of phase 
with the RE fields when it arrives at the accelerator, so that the residual electron beam 
deposits most of its kinetic energy into the accelerator’s RE fields. Then, the much lower 
energy electrons departing the accelerator are directed to the beam dump. 

8. Beam Dump 

At the end of the electron beam line, the beam dump is designed to safely and 
efficiently absorb the residual kinetic energy of the electron beam while dissipating the 
heat and shielding the radiation that is generated by collisions between the electrons and 
the beam dump material. The complexity of the beam dump is proportional to the 
remaining power in the electron beam. For a “low energy” beam, the beam dump might 
just be a copper block or Faraday cup that is air cooled, but for a “high energy” beam, the 
beam dump becomes a complex and intricate system in its own right. 

Since the beam can often have megawatts of power remaining in high energy 
systems, the heating of the beam dump and breakdown of the dump material can be a 
significant concern. Some designs feature a long conical hole where the beam strikes to 
slowly “shave” the edges of the beam and spread the deposited heat over a wider area, 
and all designs utilize cooling to prevent failure of the beam dump material. While these 
cooling systems can be quite complex in structure and function, the primary purpose is to 
safely dissipate the thermal energy deposited by the beam, and most designs utilize a 
simple water recirculation system [14]. 
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If the electrons in the beam are energetic enough, the impact of the beam can 
generate serious radiation hazards to equipment and personnel. Obviously, the “high 
energy” electrons themselves can cause considerable biological damage, but the 
collisions involved can generate neutrons and even gamma radiation in some cases. 
Therefore, the beam dump is usually surrounded by massive amounts of shielding 
material (often concrete blocks) to contain the radiation generated during operation. In 
many cases, the beam dump is the one component that is physically isolated from other 
portions of the system, e.g., in its own room, to allow for easy personnel access to the 
systems without exposing them to the residual radiation in the dump. In shipboard 
applications, the beam will enter the dump at energies low enough so that neutrons are 
not produced. 
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B. FEL RESONANCE CONDITION 

For optimum energy exchange between the electrons and the optical wave, the 
electrons’ velocities must be such that over the distance of one undulator period an 
electron will fall behind a photon traveling at the speed of light by approximately one 
optical wavelength. As an electron passes through one undulator period, approximately 
one wavelength of light passes over it. 

This “electron-photon” race is conceptually illustrated in Figure 9 with a red 
electron, blue photon (of wavelength X), and green undulator period (ko). Due to its 
slightly lower speed, the “red” electron only travels a distance of kg, while the blue 
photon travels a total distance of ko + k. In order to determine the wavelength of light 
emitted, the velocity difference between photon and electron (c(l - Pz)) must be 
multiplied by the time it takes the electron to cover one undulator period (ko/ ). 
Assuming relativistic electrons, these substitutions and some simplifications produce the 
following relations between the optical wavelength (k), the undulator period (ko), the 
dimensionless undulator parameter (K), and the Lorentz factor (y »1) 



where 


(III) 


( 11 . 2 ) 


where e is the electron charge. Me is the electron mass, c is the speed of light, and Brms is 
the root mean squared value of the undulator magnetic field. Equation (II. 1) 
demonstrates how an FEL can be “designed” or tuned for different wavelengths, by 
changing ko, K, or y. 
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Eigure 9. The “Electron-Photon Race.” Erom [15] 
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C. FEL PENDULUM EQUATION 

The electron dynamics within a free electron laser can be properly described by the 
solutions of the Lorentz force equations that account for the magnetic field of the undulator 
and the electric and magnetic fields of the optical pulse. The dynamics are controlled by 
relativistic but non-quantum mechanical effects—only classical mechanics and 
electromagnetism are applied. An “ideal” helical undulator has a magnetic field () of 


=^«,[cos(^„z„),sm(^„z„),0] , 

(11.3) 

where Bm is the undulator magnetic field strength, ko is the undulator wavenumber 

{ko=2n/Xo), and z^ is the distance along the undulator axis. The associated optical electric 

and magnetic fields, in cgs units, are 


= E[cosip,-smip,0] , 

(11.4) 

=E’[sinV’,cosV’,0] , 

(11.5) 

where i// = kzu- cOopt + (p, E is field amplitude in cgs units, (p 

is the phase, k is the 

wavenumber {k = 2k/X), X is the wavelength, and cOop is the optical wave frequency. For 

a given electron velocity (vg), the relativistic Lorentz force equations 

are 

at m^c ' ' 

(11.6) 

7= 7= ,and 

at m^c ' ' 

(11.7) 

r" = \-W , 

(11.8) 


where (3 = . 

Assuming that the electrons are traveling at relativistic velocities and that they 
experience perfect injection into helical orbits, Pz ~ 1 and the transverse motion of the 
electrons in the undulator is 
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(II.9) 


^±=(“%)bos(^„z„).sin(i„z„),0 . 

After combining equation (II.4) with equation (II.7), it can be shown that 

pKF 

i = - -cos(C + 0 ) , (11.10) 

where C = (^ + ko)zu - (Oopt is the “electron phase”. Since k»ko, the electron phase 
dsecribes the electron’s microscopic position within an optical wavelength as it passes 
through the undulator. By defining the “electron phase velocity” as i)(x)=d^dx=f [(A: + 
ko)^z- k\, we can derive the PEL “pendulum equation” 

o oo 

z/ = C = Hcos(C + 0) , (11.11) 

where v denotes a derivative of the electron phase velocity with respect to the 
dimensionless time, r = is the dimensionless time, L is the length of the undulator, 

N is the number of undulator periods, and | a | = 4nNeKLE/ ylnieC^ is the dimensionless 
optical field amplitude. The laser has weak optical fields that do not over bunch electrons 
when I a I < 71, but when | a | » 7i, the laser has strong fields and is near saturation, 
where electron over-bunching takes place. The specific phase-space path of each 
individual electron can be identified by 

j /2 ^ul + 2 |a|[sin(C + 0) - sin(C + 0)] , ( 11 . 12 ) 

where and define the initial electron coordinates in phase space, and the separatrix 

is given by = 2|a|[l-|-sin((^-|-0)]. The amplitude of the closed orbit in phase space 

• ol iX 
is2 . 
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D, FEL WAVE EQUATION 


In the presence of an electron beam, the full wave equation is 


V 


1 d 


2 t 


c' de 


c 


(11.13) 


where is the distance from the beam axis, A(x^,t) = cosV^,—sm'i/^,0] is the 


optical vector potential, and is the current density due to the undulator’s oscillating 
magnetic fields. The laser’s electric and magnetic fields are related to the vector 
potential, A , by Maxwell’s equations 


1 dA 

-and 

c dt 

(11.14) 

VxA . 

(11.15) 


After combining Equations (11.13) through (11.15) above with Equation (II.9) and 
assuming slow variations in the field’s phase and amplitude, we derive that the field 
evolves according to 


d 

dt 



InKecp^ 


(11.16) 


where Pe is the electron density and (...) is the average of all electrons. If we then 

2 2 * 

multiply both sides of equation (11.16) by a factor of diiNeKE/ j mgC , it can be written as 



d 

dr 


a{x,t) = -[je , 


(11.17) 


where the dimensionless current density (j) is defined as 



In equation (11.17) above, we ean traee the effeets of eaeh term on the diffraetion, 
gain, and bunehing in an PEL. The term operating on the dimensionless laser field 

amplitude, a , deseribes the diffraetion of the laser beam. When diffraetion is signifieant, 
the left side of equation (11.17) ineludes a double derivative term to aeeount for the 
ehanging amplitude and phase; but when diffraetion is small, the wave equation 
simplifies to equation (11.19) where only the dimensionless time derivative term of a 
remains 


°a = -j{e-'^) . (11.19) 

The dimensionless eurrent density, j, deseribes the eoupling between the laser 
light and the eleetron beam and therefore the gain of the free eleetron laser. When the 
dimensionless eurrent is small (y < tt ) the FEE eoupling and gain is small, but when 

7 ^ TT , the eoupling and gain are large. The bunehing of the eleetron beam is deseribed 

by the final term {e^'^ ^. When the eleetrons are randomly distributed in the overall 

average and eoupling is small, so that the wave does not evolve, a ^ 0 . 

As the laser begins operation, a feedbaek loop develops between the wave 
equation (11.19) and the pendulum equation (II. 11). While traveling along their phase 

spaee paths, the eleetrons begin to buneh, this eleetron bunehing ehanges the 'j term 

in equation (11.19), and therefore the dimensionless time derivative of the laser field 
amplitude also ehanges. The laser field, amplitude or phase, ean inerease or deerease 

depending on the sign of 'j. An inereasing a ean also drive an inerease in the 

eleetron phase aeeeleration from equation (II. 11). This inerease in phase veloeity fuels 
further bunehing, ereating a feedbaek loop that will drive gain exponentially higher. This 
field growth eventually ean fall in strong optieal fields as the eleetrons move through 
their paths in phase spaee and begin to absorb energy baek from the optieal beam. 
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E. 


EEL PULSE STRUCTURE 


The MW-class PEL is often referred to as a high average power laser, while most 
other lasers that produce the same average power are referred to here as high-peak power 
lasers. An important difference between the FEE and other lasers of similar power is its 
pulse stmcture. High-peak power lasers can produce peak powers in the TW range, but they 
maintain a low pulse repetition frequency (PRF) in the tens of Hz. High average power 
lasers, like a MW-class FEE, produce a peak power in the GW range, but their PRF is almost 
1 GHz. Figure 10 displays the micro-pulse structure for a MW-class FEE with repeated 
micro-pulses of over 1 GW at approximately 700 MHz PRF producing an average power of 
about 1 MW. Eater analysis, in Chapter VII, will utilize the ratio of peak-to-average power 
in a comparison of laser and solar intensities. It is important to realize that the FEE peak 
power is over 1000 times greater than the average power, and that the PRF is highly precise. 



TIME 

(Not to Scale) 

Figure 10. MW-Class FEE Macro-Pulse Structure 
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III. LASER PROPAGATION 


The application of a high-energy, ground-based laser for power beaming, 
extraterrestrial illumination, or orbital debris clearance requires an understanding of beam 
propagation. Transmitting a high power beam through a possibly thick and fluid 
atmosphere over extreme ranges can affect the beam quality and lead to unintended 
effects at the target. Random thermal and density gradients within the air cause random 
fluctuations in the refractive index of the atmosphere that degrade the propagation of an 
optical wave. These effects cause the twinkling of stars that we see at night and limit the 
resolution ability of earthbound telescopes to a few seconds of arc, but they can also 
degrade the spatial coherence of a laser beam as it propagates through the atmosphere 
[16, p. 1]. The loss of spatial coherence will ultimately limit the ability to focus the laser 
beam, and it can drastically lower the power on target of laser systems. In most cases, 
these losses can be mitigated through the use of an adaptive optics system for the 
director. Also, like astronomical observatories, ground-based lasers utilized in space 
applications could avoid most of the detrimental effects of atmospheric propagation by 
operating at high altitude, but that would require the construction of an associated 
electrical generator. 

The fundamental propagation effects faced by a ground based laser are: 
diffraction, scattering, absorption, turbulence, and thermal blooming. With the notable 
exception of thermal blooming, all of these are linear effects, and the relative importance 
of each loss mode varies with application. As described by Andrews, some common 
types of space and laser applications are cited below with a brief description of their 
primary atmospheric effects after [16, p. 11-13]: 

• Satellite-ground: Laser communications to the ground from a satellite are 
disrupted by the atmospheric turbulence near the ground, but for most of 
the path the beam passes through little turbulence. Because the 
propagation distances are so long, the beam is very wide by the time it 
encounters the atmospheric layer. Hence the beam is mostly disturbed by 
spatial phase fluctuations, thereby limiting heterodyne efficiency in 
coherent detection. The primary concerns for downlink propagation paths 
are scintillations and angle-of-arrival fluctuations. 
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• Ground-satellite: A transmitted laser beam for eommunieations from the 
ground to a satellite is disrupted by atmospherie turbulenee near the 
ground and, thus, near the transmitter. Beeause most of the propagation 
path lies beyond the atmospherie layer, there is a long propagation path in 
free spaee that is dominated by free-spaee diffraetion. The primary 
coneerns for an uplink path are scintillations, beam wander, and beam 
pointing. 

• Aircraft-satellite and satellite-aircraft: These two communication paths 
are similar to the ground-satellite and satellite-ground paths described 
above. Although the aircraft is usually above much of the natural 
atmospheric ground-induced turbulence, aircraft boundary layer effects 
due to platform speed need to be addressed. 

• Imaging: Problems associated with imaging through the atmosphere are 
similar to those associated with propagation. For example, the “dancing” 
of an image in the focal plane of an imaging system is mathematically 
equivalent to the wander of a beam focused at the object by the same 
optical system. The resolution of a long-exposure image is equivalent to 
the long-term beam spread of a focused beam. The short-exposure 
resolution is equivalent to an optical heterodyne receiver that employs tilt- 
correction of the signal or of the local oscillator. Among others, adaptive 
optics systems are widely used today to provide turbulence-compensation 
techniques to improve image quality. 

A, DIFFRACTION 

Diffraction effects are a consequence of the wave nature of light that cause beam 
spreading as the wave propagates away from the beam waist, i.e., that part of the beam 
where the beam has the smallest diameter. This results in decreased beam intensity and 
constantly changes the phase front radius of curvature of the propagating wave. 
Diffraction is the only mechanism that is not linked to the presence of matter and, thus, in 
the approximate vacuum of near-Earth and interstellar space, it is responsible for the 
primary loss in intensity of the optical beam. Inside the atmosphere, the loss of beam 
intensity due to diffraction can be magnified significantly by atmospheric turbulence. 

When evaluating diffraction effects, there are two distinctive models available. In 
Fraunhofer, or far-field, diffraction, the point of observation is far enough from the beam 
waist so that waves arriving at the point of observation may be considered plane waves. 
When this is not the case and the curvature of the wavefront must be taken into account. 
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Fresnel, or near-field, diffraction must be applied. In the far-field approximation, the 
diffraction pattern (or beam spot size) changes uniformly in size with distance from the 
beam waist. In the near-field, the diffraction pattern is more complex [17, p. 324]. 

In Figure 11, the growth in cross sectional area or spot size of a Gaussian beam is 
graphically depicted showing the radius at the beam waist (wo), Rayleigh range (zr), 
beam divergence angle (0), range from beam waist to the target (z), and beam radius at 
the target (w). In the fundamental mode of a Gaussian beam, the intensity of the beam 
has a normal, i.e., Gaussian, distribution. In this case, the beam spot size is simply the 
area of the fundamental mode at a specific point along the beam’s axis of travel. 




Figure 11. Collimation of a Gaussian Beam. After [ 18] 


Assuming a Gaussian beam, the change in beam spot size due to diffraction in the 
far field can then be approximated based on the wavelength of the optical wave (A,), the 
range (z), and the beam waist (wo) by the equation 


W'(z) = K 



(III.1) 


For large z, the beam waist can be approximated by 
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When discussing the ability of a laser to retain a focused beam pattern, scientists 
utilize a parameter called the Rayleigh range (zr). The Rayleigh range is defined as the 
distance from the beam waist to the transverse plane where the beam area has doubled. It 
is the characteristic distance over which a beam area expands due to diffraction. 
Similarly, the beam divergence angle (0b) is used to determine the beam spot size at the 
target. The formulas for the Rayleigh range and the beam divergence angle are 


_ 7TW^ 

R A and 


(IIL3) 


0b = 


X 




(IIL4) 


The beam’s spot size (As) at a distance z from the waist is 


2 2 


4 = 


A z 


TTW^ 


2 = and 


(III.5) 


the intensity (7) is given by 


1 = 


P 


(IIL6) 


B, SCATTERING 

When any electromagnetic wave travels through a medium, it can interact with 
scattering centers in the medium, such as water vapor, suspended aerosols, or other 
particulates. These scattering centers can redirect or disperse the energy of the 
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electromagnetic wave, and in the specific case of a laser beam, scattering results in the 
loss of intensity at the target. Even on a clear day, scattering can reduce transmission in 
the visible spectrum to 20% for a 16 km path [19, p. 24], 

Mie scattering is a complex theory that takes into account the size, shape, 
refractive index, and absorptivity of the scattering particles. In the discussions below, 
Mie scattering is a term that is applied to the scattering process when the scatterer is 
comparable in size to the wavelength of the incident radiation. Rayleigh scattering and 
non-selective scattering are special cases of Mie scattering. Rayleigh scattering is the 
scattering of electromagnetic radiation by particles much smaller than the wavelength of 
the radiation; i.e., individual atoms or molecules. Nonselective, or geometrical, 
scattering occurs when the scattering center is much larger than the wavelength of the 
incident light, such as with haze, fog, or rain. Comparing the size of the scattering center 
and the wavelength of the incident radiation, as shown in Table 1, allows for 
simplifications of Mie theory into Rayleigh scattering, when the scatterer is small, and 
non-selective scattering, when the scatterer is large. 


Type of Scattering 

Size of Scatterer 

Rayleigh Scattering 

Larger than an electron 

but smaller than X 

Mie Scattering 

Comparable in size to X 

Non-selective Scattering 

Much larger than X 


Table 1. Types of Atmospheric Scattering. From [19, p. 24] 


For all types of scattering, the bound electrons experience a displacement due to 
the harmonic oscillation of the incident electric field. The response of the electron to this 
driving force is a function of the driving frequency and the natural or resonant frequency 
of the oscillator [17, p. 303]. The oscillating dipoles re-radiate or scatter their energy in 
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all directions except along the dipole axis itself, but each form of scattering results in a 
different scattering pattern. As the scattering center grows in size with respect to the 
wavelength of radiation, the energy reemitted by the dipole begins to emit preferentially 
in the direction of travel of the incident light. This is graphically illustrated in Figure 12. 


Rayleigh Scattering 



Mie Scattering 



Direction of incident light 


Mie Scattering, 
larger particles 



Figure 12. Rayleigh, Mie, and Nonselective Scattering. After [20] 


1, Rayleigh Scattering 


Rayleigh scattering is the scattering of electromagnetic radiation by particles 
much smaller than the wavelength of the radiation; i.e., individual atoms or molecules. In 
Rayleigh scattering, the distance between scattering centers allows them to act 
incoherently, and therefore, their net irradiance is a sum of their individual irradiances. 
From [21, p. 448], the radiated power can be shown to be proportional to the fourth 
power of the frequency or inversely proportional to the fourth power of the wavelength of 
the incident radiation. This is the Rayleigh Scattering Law 



4 2 

^^o^fPo 

IIttc 


(IIL7) 


where Pdp is the power radiated by the dipole, po is the permeability of free space, oof is 
the frequency of the incident radiation, and po is the maximum value of the atom’s dipole 
moment. Therefore, as shown in Figure 13, oscillating dipoles scatter more energy in the 
high frequency, short wavelength region of the optical spectrum than in the low 
frequency, long wavelength region. 
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WAVELENGTH (pn) 

Figure 13. Rayleigh Scattering. From [19] 


The sun emits radiation across most of the electromagnetic spectrum, but its 
output is not uniform at all wavelengths (covered in more detail in Chapter V). Within 
the optical spectrum, the peak of the sun’s output is in the yellow wavelengths, and when 
viewed directly, the sunlight appears to be yellowish-white in color. Inside the Earth’s 
atmosphere, sunlight experiences Rayleigh scattering due to the nitrogen and oxygen 
molecules in the air. If there were no atmosphere, the sky would appear black, but inside 
the atmosphere, the short wavelength, blue light is scattered randomly and arrives at the 
observer from all directions, making the sky appear blue. Viewing the sun at its zenith, 
the sunlight travels through a relatively thin portion of the atmosphere and only a fraction 
of the light is scattered, making the sky appear blue while the sun appears yellow or 
white. Viewing the sun low on the horizon, the sunlight travels through a relatively long 
and dense section of the atmosphere and the sun appears red or orange. The preferential 
scattering of shorter wavelengths in the daytime sky is a clearly visible demonstration of 
the 1 /effect in Rayleigh scattering. In the design of a ground based laser system for 
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space applications, the implication is clear that, within the optical spectrum, longer 
wavelengths will experience less loss due to seattering [17]. 

2, Mie Scattering 

Mie scattering occurs when the scattering eenters are comparable to the size of the 
incident wavelengths, such as with aerosols or small water droplets. In Rayleigh 
scattering, the spatial variation of the electric field over the molecular charge distribution 
could be ignored due to the relative sizes of the wavelength and molecule. When the 
radius of the molecule is approximately equal to the wavelength of the incident radiation, 
the variations of the electric field over the molecule become significant, and they can 
drastically affect the directionality of the dipole radiation. Despite aerosol particle 
concentrations that are orders of magnitude lower than molecular concentrations in the 
atmosphere, Mie scattering in the optical wavelengths far exeeeds the attenuation due to 
both Rayleigh seattering and ozone absorption [19, p. 29]. As shown in Figure 14, the 
Rayleigh scattering coefficient only beeomes eomparable to the total aerosol extinction, 
i.e., aerosol scattering plus absorption, for wavelengths shorter than about 0.3pm, 
whereas aerosol seattering and absorption dominate at longer wavelengths. 



Figure 14. Rayleigh and Aerosol Attenuation. From [22] 
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3. 


Nonselective Scattering 


Nonselective, or geometrical, scattering occurs when the scattering center is much 
larger than the wavelength of the incident light, such as with haze, fog, or rain. Haze 
refers to the small particulates of salt, dust, and combustion products suspended in the air 
with radii less than about 0.5 pm. In conditions of high humidity, these small particles 
serve as nucleation sites for the condensation of water molecules, and this condensation 
can rapidly increase the overall surface area of the particulate. When these particulates 
increase in size to exceed 0.5 pm, the haze becomes fog if touching the ground, or clouds 
if suspended above it. 

As the particles continue to grow, they become raindrops, which are many times 
larger than the wavelength of a laser beam. In the presence of these large scattering 
centers, the wavelength dependence on scattering is eliminated, and the scattering 
coefficient varies purely as a function of the drop size [19, p. 41]. Thus, all wavelengths 
scatter equally from clouds, making them appear white. In a heavy rain, laser beam 
transmission through the atmosphere can be dramatically affected. Rainfall rates for four 
different rain conditions and the corresponding transmittance, due to scattering only, of a 
1.8 km path are shown in Table 2. 


Condition 

Rainfall (cm/hr) 

Transmittance 

of a 1.8 km path 

Light Rain 

0.25 

0.88 

Medium Rain 

1.25 

0.74 

Heavy Rain 

2.5 

0.65 

Cloudburst 

10 

0.38 


Table 2. Transmittance of a 1.8 km Path Through Rainfall. From [19, p. 43] 
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In Figure 15, the attenuation eoeffieients for a typieal dense fog and light rainfall 
is plotted versus wavelength from 0.6 to 15 pm. To allow for eomparison, the absorption 
eoeffieients of a layer of water eontaining the same amount of liquid water as the rain are 
also presented. For all wavelengths shown, the attenuation of rain is relatively eonstant 
and mueh less than that of a dense fog. It is immediately apparent that any ground based 
laser system intended for spaee applieations should operate at a high altitude and arid 
ground site to minimize the losses due to seattering regardless of the wavelength in use. 



0.6 1 1.5 3 6 10 16 


WAVELENGTH (jiJii) 


Figure 15. Attenuation of IR Radiation due to Rain and Fog. From [19, p. 44] 


C. ABSORPTION 


Absorption is the proeess by whieh the energy in an eleetromagnetie wave is 


attenuated in passing through a medium by the eonversion of the energy into heat or other 
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forms of molecular energy. The capacity of a specific medium to absorb radiation is 
dependent on several factors such as the wavelength of the radiation, the thickness of the 
medium, the atomic and molecular properties of the medium, and the temperature and 
density of the medium. A perfectly transparent medium would permit the passage of a 
laser beam without any loss of intensity, while a perfectly opaque medium would 
completely block the beam. In reality, there is no material that is perfectly transparent or 
opaque across the entire electromagnetic spectrum. Thus, it is vitally important to 
consider the absorption spectrum of the intervening medium when attempting to 
propagate an electromagnetic wave through the medium. Most importantly, the amount 
of absorption determines the likelihood of thermal blooming. This process is nonlinear, 
and it can limit the power transmitted over long distances through the atmosphere. 
Specifically, any ground based laser must find a “window” of atmospheric transparency 
at some appropriate wavelength in order to maximize the transmission of beam power 
through the atmosphere [23]. 

Earth’s atmosphere is composed of 78% nitrogen, 21% oxygen, and 1% of argon and 
trace elements. Nitrogen (N 2 ) and oxygen (O 2 ) have no dipole moment and therefore do 
not exhibit molecular absorption bands. In the other gases, water (H 2 O), carbon dioxide 
(CO 2 ), and ozone (O 3 ) are the most significant contributors to atmospheric absorption. 
Water vapor varies widely, but averages about 0.5% of the atmosphere. Ozone is found 
in relatively large concentrations at high altitudes near 30 km, but the concentration of 
ozone near sea level is negligible [19]. 

Figure 16 is a chart of the atmospheric transmittance, including scattering and 
absorption effects, measured over an 1820 m horizontal path at sea level. Table 3 lists 
the wavelength intervals where transmittance is relatively high. For efficient 
transmission, the laser wavelength should be selected in the center of these windows as 
the transmission is decreased at the window boundaries. Figure 16 and Table 3 are useful 
in examining the general effects of absorption on atmospheric transmission over near¬ 
horizontal paths, but any conclusions drawn from them would not necessarily apply for a 
ground based laser system transmitting vertically through the atmosphere. The FEE 
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bandwidth is estimated at 0.1 %, or about 10'^ |a.m for an PEL operating at 1 pm, and it 
would be expected to “fit” into any of these windows [15]. 



Figure 16. Atmospheric Transmittance. From [19, p. 15] 


Window Number 

Window Boundaries (pm) 

Low High 

I 

0.72 

0.94 

II 

0.94 

1.13 

III 

1.13 

1.38 

IV 

1.38 

1.90 

V 

1.90 

2.70 

VI 

2.70 

4.30 

VII 

4.30 

6.0 

VIII 

6.0 

15.0 


Table 3. Wavelength Regions of Atmospheric Windows. From [19, p. 16] 
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D, TURBULENCE 


Treating the atmosphere as a viseous fluid, the motion of the air ean be deseribed 
by fluid meehanies as either laminar or turbulent. In laminar flow, the veloeity is 
uniform, or ehanges in a regular pattern, and mixing does not oeeur. In turbulent flow, 
dynamie mixing of the air leads to the ereation of turbulent eddies within the air. Like a 
layered onion, the Earth’s atmosphere is relatively dense within the flrst 20 km, but 
extends in varying layers to approximately 700 km. In a region within about 2 km of the 
Earth’s surfaee, ealled the atmospherie boundary layer, atmospherie dynamies are 
dominated by the viseosity of the air and the heat exehange interaetion with the Earth’s 
surfaee, and the atmosphere ean be turbulent [16]. 

Atmospherie turbulenee is a linear effeet that oecurs when air paekets at differing 
temperatures are mixed by wind and eonveetion. Over time, the air will reaeh a thermal 
equilibrium as the turbulenee eells break down into smaller eddies. In the interim, the 
turbulent mixing produees random density fluetuations, and therefore, ehanges in the 
index of refraetion of the air. 


This turbulent flow of air, temperature fluetuations, and the assoeiated density 
fluetuations within the propagation path of a laser beam ean eause variations in the index 
of refraetion along the path length. These small and loealized variations in the index of 
refraetion ean then eause phase fluetuations within the laser beam. The phase 
fluetuations aet as random lenses, and they result in amplitude fluctuations across the 
beam wave front. Wave front distortions caused by atmospheric turbulence can lead to 
the broadening of a laser beam, random changes in the beam direction (“beam wander”), 
and intensity fluctuations within the beam profile (“scintillation”). 

The defining parameter to describe the beam disturbances caused by turbulence is 
the refractive index structure coefficient (Cn). It is a function of the pressure and 
temperature difference between two points separated by distance r, and it is given by 



79x10 



(IIL8) 
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where the temperature strueture parameter (Cj) is 


C^ = 

T is the average air temperature in degrees Kelvin, Ti and T 2 are the temperatures at the 
points of interest in degrees Kelvin, Pp is the atmospherie pressure in millibars, and r is 
the distanee in em between points. Measurements of €„ inside the atmospherie boundary 
layer have demonstrated that it reaehes a minimum about one or two hours before sunrise 
and after sunset [19, p. 49]. Typieal values of the refraetive index strueture are shown in 
Table 4. 




(III.9) 


Strength of Turbulence 

C„ 

Strong 

5 X lO'^ 

Intermediate 

4 X 10'^ 

Weak 

8 X lO'*^ 


Table 4. Typieal Values of C„. From [19] 


1. Beam Wander 

Beam wander deseribes the random motion of the beam’s spot at the target due to 
turbulenee effeets that are generally larger in size than the beam diameter. In general, 
beam wander has been demonstrated to be independent of the beam wavelength, and 
elosely follows the value of C„ [19, p. 54]. It beeomes a major problem for applieations 
at long range when the laser must stay fixed on a relatively small target. Thus, for a 
ground based laser beaming power out of the atmosphere, a fast and highly aeeurate 
optieal traeking system will probably be required to keep suffieient intensity on target. 


40 





2 . 


Scintillation 


Scintillation describes the “twinkling” of stars or the “flieker” of a laser beam as it 
arrives at the target. When the air flow moves small index-of-refraetion inhomogeneities 
through the path of the beam, the dismption of the beam wavefront ean lead to intensity 
variations within the beam profile and result in “hot spots” at the target. Although the time 
averaged intensity on target remains fairly eonstant, seintillation ean be harmful in laser 
applieations. Loeal “hot spots” within the beam ean lead to thermal blooming and in power 
beaming applieations the higher intensities may exeeed the solar eells’ power eonversion rate. 

3, Turbulence Induced Beam Spreading 

In addition to diffraction effects, a laser beam propagating through the turbulent 
atmosphere experiences a broadening due to the loss of spatial coherence of the wave. 
Sometimes referred to as “beam breathing,” this effect is similar to diffraction, and it can 
increase the beam spot size at the target. In practice, observed spot sizes can be twice as 
large as those predicted by diffraction theory alone [19, p. 61]. 

E. THERMAL BLOOMING 

For high power beams, the molecular absorption of the beam’s energy by the 
atmosphere can lead to temperature, density, and index-of-refraction changes within the air. As 
the beam continues to travel through air, it can experience a nonlinear, defocusing or 
“blooming” effect. 

Figure 17 shows the evolution of a Gaussian beam and its temperature (T), density (p), 
and index of refraction (n) as the beam suffers the effects of thermal blooming. No turbulence 
is included in this example. The lower half of Figure 17 shows the initial irradiance, 
temperature (To), density (po), and index of refraction (no) profile of a Gaussian beam. At the 
center of the Gaussian beam, where the intensity is the highest, the absorption of laser energy 
and subsequent release of heat results in a radial temperature profile with the hotter, less dense 
air in the center of the beam and the colder, more dense air at the edges of the beam. These 
density changes lead to a change in the index of refraction which tends to deflect the power of 
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the beam radially away from the beam center. Thus, a Gaussian beam at transmission arrives at 
the target in a ring pattern, as shown at the top of Figure 17. 
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Figure 17. Thermal Blooming of a Gaussian Wave. From [19, p. 69] 
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Figure 18 is an illustration of a laser beam suffering from thermal blooming. In 
this figure, a laser beam of power (P) and cross sectional area (A) propagates through an 
air column of total length (R). During its propagation, the laser beam passes through a 
stagnation range (Rg) representing a column of stationary air. It is within this stagnant 
column of air that thermal blooming occurs and begins to spread the beam. An 
estimation of the temperature increase per second is given by 


Ar aRP 
_ ± _ 

At pCpAR^ ’ 


(III. 10) 


where AT is the change in temperature given in K, At is the change in time given in 
seconds, a is the absorptivity of air (~ 0.0041 km'* at 1pm [25]), Rg is the stagnation 
range, P is the power of the laser, p is the density of air (~ 1.2 kg/m^ at sea level), Cp is 
the specific heat of air (~1 kJ/kg-K), and A is the cross-sectional area of the beam. [24] 
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Figure 18. Thermal Blooming Graphic. From [24] 

The increase in temperature leads to a change in the index of refraction of the air 
according to the formula 

This small change in the index of refraction ( An = — ) results in a small velocity 

difference Av for light in the middle part of the beam, which is warmer, and the edge of 
the beam, which is cooler. Thermal blooming is created when this velocity difference 
allows the wave front to curve by as much as T over the path length Rg. So, the index 
change An causes thermal blooming when 

An^-y^^ , (III. 12) 
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where An is the dimensionless index of refraetion and X is the laser wavelength in meters. 
By eombining Equations (III.l 1) and (III. 12), we ean obtain an equation for the ehange in 
temperature that would lead to thermal blooming 


AE: 


lO^’A 

R 


in degrees K. 


(III. 13) 


Combining Equations (III. 10) and (III. 13) results in an equation for the time 
interval until thermal blooming begins 




(2xlOV/m^)A^ 

aRf 


(III. 14) 


where AIb is the time until thermal blooming begins in seeonds. [24] 

Equation (III. 10) approximates the temperature ehange per seeond due to the 
absorption of beam power, and it ean be applied to the MW-elass EEE. Assuming a 
power of 1 MW, beam area of about 3.14 m , and a stagnation range of 10 km, it prediets 
a value of 1.08 x 10'^ K/see. This is eneouraging, until Equations (III.13) and (III.14) 
show that thermal blooming might oeeur for a temperature ehange as small as 10'"^ K, and 
it eould begin as quiekly as 0.09 seeonds after initial illumination. Obviously, this would 
not allow for sufficient propagation of the beam over long ranges. 

Appearing in the numerator and denominator of Equation (III. 10), the stagnation 
range does not affect the rate of temperature change due to absorption. However, it can 
have a large effect on Equations (III. 13) and (III. 14). In the case presented above, the 
initial estimation that the entire atmosphere is stagnant is extremely conservative. There 
are many variables that affect thermal blooming including the atmospheric density, the 
velocity of the air, the range of propagation, and the angular rate of the beam as it tracks a 
target. Eor the power beaming applications discussed in Chapter V, the laser would 
sweep across the sky as it tracks a target at an angular rate of approximately 1/2 degree 
per second. Theoretically, the beam might only experience thermal blooming over very 
short stagnation ranges close to the director mirror. If this stagnation range was about 
100 m, then thermal blooming might not begin until approximately 9 seconds after the 
first illumination. 
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Likewise, some of the effects of thermal blooming could be mitigated by locating 
the ground-based laser in a dry climate at high altitude to reduce the atmospheric 
absorption, propagation range, and stagnation range of the laser beam. However, for all 
of the topics within this thesis, thermal blooming effects are mostly ignored to obtain a 
first approximation of the feasibility of these applications. While the complex issues 
involved in the atmospheric propagation of the laser are not trivial, over thirty years of 
research within the military indicates that a high-powered laser can be successfully 
propagated through the atmosphere. The military has been interested in FELs since 1978, 
when “the Defense Advanced Research Projects Agency (DARPA) concluded that no 
other high-power laser could achieve the optical beam quality necessary to focus the 
beam on a distant (thousands of km) target [26, p. 2-1].” These topics obviously require 
further study, but for the purposes of this thesis they will be generally excluded from 
consideration. 
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IV. ORBITAL MOTION 


The application of a ground based laser to targets outside of Earth’s atmosphere 
requires a detailed understanding of the variables involved. Both the beam’s propagation 
through the atmosphere and the motion of the target are critical to the design of an 
integrated system. In most of the applications discussed here, the target will be a satellite 
or satellite debris in Earth orbit; therefore, an understanding of orbital mechanics and the 
space environment will be required. This chapter will describe the basic laws of orbital 
motion and describe the forces that act on an orbiting object. A general comprehension 
of the relative magnitudes of the perturbing accelerations is conceptually important. 
Chapter VI will attempt to utilize some of the perturbations generated or altered by a 
ground-based laser to modify the orbits of space debris. Specifically, the laser may be 
able to increase the aerodynamic drag forces on orbital debris by altering the orbit of the 
debris through vaporization of a small amount of surface material, break-up of the debris 
into smaller particles, and radiation pressure. While most of the other perturbation 
effects will not be specifically calculated in later work, many of them are utilized within 
the code of the Satellite Tool Kit (STK) program. 

A, BASIC THEORY 

Despite their limited tools, ancient astronomers accomplished many insightful 
discoveries regarding the motion of celestial bodies, but prior to the 1600s, most of these 
discoveries amounted to cataloging and describing specific interactions, such as 
predicting the changes in lunar phases, solar eclipse cycles, and the phases of Venus and 
other nearby planets. Building on the Copernican, heliocentric model of the solar system 
and Tycho Brahe’s very precise observational data, Johann Kepler produced a work, 
Astronomica Nova (New Astronomy), that would finally capture the kinematics of 
planetary motion and lay the foundation for the field of orbital mechanics. Any review of 
two-body, Keplerian, mechanics requires a familiarity with these laws. 
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Kepler’s Laws: 

1) Law of Ellipses; An orbit is an ellipse, or eonie seetion, with a eentral 
body at one foeus. 

2) Areal Law: An orbiting body’s radius veetor from a eentral body sweeps 
out equal areas in equal times. 

3) Harmonie Law: The square of an orbiting body’s revolutionary period is 
proportional to the eube of the satellite’s mean distanee from the eentral 
body. 

In the early 1700s, Galileo Galilei provided an approximate value for the 
aeeeleration due to gravity, stating that two bodies dropped simultaneously would fall at 
the same veloeity independent of their masses. While Kepler deseribed kinematie motion 
that details the motion of objeets without eonsideration of the eauses leading to the 
motion, Galileo took the first steps into orbital dynamies and attempted to deseribe the 
eauses of orbital motion. Synthesizing Kepler’s laws and Galileo’s terrestrial meehanies. 
Sir Isaae Newton diseovered the universal law of attraetion and proposed the three 
fundamental laws of elassieal meehanies. When eombined with Kepler’s laws, Newton’s 
seeond law and his universal law of gravitation are the starting points for any modern 
study of orbital motion [27, p. 29-31]. 

Newton’s Laws: 

1) The Law of Inertia: Every body eontinues in its state of rest, or of uniform 
motion in a straight line, unless it is eompelled to ehange that state by 
forees impressed upon it. 

2) The Lundamental Law of Dynamies: The ehange of motion is proportional 
to the motive foree impressed and is made in the direetion of the right line 
in whieh that foree is impressed 

3) The Law of Aetion and Reaetion: To every aetion there is always opposed 
an equal reaetion; or, the mutual aetions of two bodies upon eaeh other are 
always equal and direeted to eontrary parts. 
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B, 


GRAVITATIONAL THEORY 


Despite the intrieate models that have been developed to eompute the motion of 
satellites to the required aeeuraeies for various applieations, the main features of their 
orbits may still be deseribed by a few simple and reasonable approximations. In the same 
way that the Sun primarily governs the motion of the planets, Earth’s gravitational foree 
exeeeds all other forees aeting on a satellite by several orders of magnitude. As a first 
approximation, Kepler’s laws of planetary motion, Newton’s classieal meehanies, and 
Newton’s law of gravity may be applied to a satellite’s orbit about the Earth. Ignoring, 
for the moment, the perturbative forees upon a satellite due to their relative 
insignifieanee, the “two-body” solution to orbital motion ean be developed. In this ease, 
the satellite’s mass is negligible when eompared to the mass of the Earth, and the Earth 
and satellite are assumed to be spherieally symmetrie and of uniform density. In its 
simplest form, the aceeleration (r ) of the satellite is given by Newton’s law of gravity 


r = 



(IV.l) 


3 2 

where p is Earth’s gravitational parameter (-398,600 km /s ), r is the distanee from the 
center of the Earth to the center of gravity of the satellite in km, and f is a unit vector 
along that path. [28, p. 15] 

In this case, the satellite is confined to orbit in a fixed plane, since the 
gravitational force exerted on the satellite always points towards the center of the Earth, 
hollowing the derivations outlined in [28, p. 15-19], the velocity {vsat) and period {Tsat) of 
a satellite in this simplified orbit is given by 


V 


sat 



and 




sat 



where as is the orbit’s semi-major axis, given by 


(IV.2) 

(IV.3) 
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a 


(IV.4) 



where h is the speeifie angular momentum and e is the eceentrieity of the orbit. For a 
eireular orbit where r=a. Equations (IV.2) and (IV.3) simplify to 



(IV.5) 


(IV.6) 

The equations presented above outline the aeeeleration and related variables 
eaused by a simplified, eentral potential. A more aeeurate solution of the aeeeleration 
due to Earth’s gravity ean be deseribed by treating the Earth as an oblate spheroid and 
expanding the geopotential to higher degrees of aeeuraey. 

Envisioning a summation of the gravitational attraetion between all of the diserete 
point masses that make up the Earth and the satellite’s mass, an expression for the 
geopotential emerges that has no known elosed-form solution. An expansion of the terms 
ean be aceomplished using a series of Eegendre polynomials and spherieal harmonies. 
Spherieal harmonies represent periodie boundary solutions on a unit sphere, and in this 
ease. Zonal, Seetoral, and Tesseral Harmonies attempt to deseribe the eomplex 
distribution of mass within the Earth. As shown in Eigure 19, Zonal harmonie terms are 
axially symmetrie about the earth’s polar axis of rotation, and they represent the 
flattening of the Earth due to its rotation and other latitudinal mass distributions. Seetoral 
harmonies divide the Earth along longitudinal meridians, and they take into aeeount the 
extra mass distributions in longitudinal regions. Tesseral harmonies “eheeker” the Earth 
in an attempt to model speeifie loeal regions that depart from the unit sphere [29, p. 509- 
520]. 
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Zonal Harmonics 


Sectoral Hannonics 


Teeseral Harmonice 



Figure 19. Spherical Harmonics. After [30] 


Ignoring the tidal variations in the terrestrial mass distribution and the variation in 
the direction of the Earth’s axis of rotation, the central potential can be written as 


U[r,K,^) 


oo 

/=! 

oo I 

+5] 5] (e) [Cim COS (mA) + sin (mA)] (sin$) 


/=1 m=\ 


, (IV.7) 


where r, A, and 0 are spherical coordinates related to the x, y, and z Cartesian 
coordinates by 

X = rsinAcos'h , (IV.8) 

y = rsin Asind?, (IV.9) 

z = rcosA, (IV.IO) 


R is the equatorial radius of the Earth, / and m are counting variables, the PfsinO^ terms 
are Eegendre polynomials, and J, C, and S are harmonic coefficients. In Equation (IV.7), 
the first term represented by the number one, describes only the dominant central 
potential causing Keplerian motion. The second term, with J/ and Pi, displays the 
contribution of the zonal harmonics. The third term, with C/m, Sim, and Pim, provides the 
contribution of the sectorial and tesseral harmonics. This expansion is discussed here to 
familiarize the reader with the even terms of J that can be some of the larger perturbation 
forces affecting an object at low altitude. The additional terms not discussed above, such 
as the Cii, Su, and J? terms, are negligible when modeling the actual mass distribution of 

the Earth, but the interested reader can find a full derivation and discussion in [27, p.67]. 
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Neglecting the expansion of the geopotential and the perturbing forces described 
in Section D, the fundamental equations of Keplerian motion, supported by Newton’s 
second law, can generally describe the motion of any orbiting object over a limited 
number of orbits. As discussed in Section E, the acceleration created by this basic 
Newtonian acceleration is about three orders of magnitude greater than the J 2 term in the 
expansion of the geopotential or the highest feasible perturbation due to atmospheric 
drag. As orbital passes build, the detailed geopotential expansion and perturbation effects 
must be used to account for variations in the orbital motion from that predicted by Kepler 
and Newton. 

D. ORBIT TYPES 

The orbital profile of any satellite has profound implications for its mission 
performance, service lifetime, electrical power, overall weight, thermal gradients, launch 
costs, and a variety of other critical factors. In addition to the complexities of systems 
engineering, there is a nearly infinite number of possible orbits for a satellite designer to 
consider during the design process. Among the thousands of satellites in orbit, however, 
there emerges a general catalog of useful orbits. There are four basic orbit types 
described below. 

The radiation environment near the Earth can be highly damaging to a satellite’s 
solar cells. It consists of electrons and protons trapped in the geomagnetic field, radiation 
generated by solar flares, and cosmic radiation. Within Earth’s magnetosphere, trapped 
electrons and protons within the Van Allen belts cause most of the damage to solar cells. 
Outside of the magnetosphere, the solar flare protons are the predominant cause of solar 
cell degradation. Eigure 20 shows a cross section of the Earth’s magnetosphere, 
illustrating the trapped electron and proton flux. Due to its highly damaging effects on 
solar cells, the center of the proton belt, between 6,000 and 12,000 km, is generally 
avoided by all solar powered spacecraft, creating an artificial gap between useful Eow 
Earth Orbits (EEO) and Medium Earth Orbits (MEO) [31, p.327]. 
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Figure 20. Earth’s Van Allen Belts. From [32] 


1, Low Earth Orbit (LEO) 

Any orbit with an altitude less than 1000 km is eategorized as a Low Earth Orbit 
(LEO). Due to their low altitude, satellites in LEO are subjeeted to signifieant 
atmospherie effeets such as orbit degradation due to atmospheric drag and material 
degradation due to erosion or chemical interactions with the upper atmosphere, but they 
can often avoid the high electron and proton flux of Earth’s Van Allen Belts [33]. A 
typical LEO satellite has an orbital period of about an hour and a half to a few hours, but 
it can only be seen from a particular terrestrial location for a few minutes as it passes 
overhead. For this reason, it is common to operate “constellations” of multiple LEO 
satellites spaced in appropriate orbits to provide frequent or continuous coverage of 
critical locations [34]. 

LEO offers many advantages over other orbits, but all LEO satellites benefit from 
the relatively low launch costs when compared to higher altitude orbits. With less fuel 
required to reach LEO altitudes, high orbital inclinations and even retrograde orbits are 
much more accessible for LEO spacecraft. While maintaining a higher viewpoint and 
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wider field of view than reeonnaissanee aireraft, remote sensing satellites in LEO are able 
to quiekly eover large portions of the Earth’s surfaee and still maintain better resolution 
than satellites at higher altitudes. Communieations satellites, sueh as Iridium or Global 
Star, benefit from the lower altitudes with inereased signal strength and shorter 
transmission delays. 

2. Medium Earth Orbit (MEO) 

Any orbit with an eoeentrieity less than 0.5 and an orbital altitude between 1000 
km and about 35,786 km is generally eategorized as a Medium Earth Orbit (MEO). 
MEO satellites ean have orbital periods between 2 and 24 hours. MEO is highly utilized 
by spaee-based, navigation satellites sueh as the Global Positioning System (GPS), 
Glonass, or Galileo systems, but some eommunieations satellites that eover the North and 
South Pole are also in MEO [34]. 

At higher altitudes than EEO, MEO satellites ean eover more of the Earth’s 
surfaee with fewer satellites while avoiding the atmospherie drag and deterioration of 
materials. Similarly, MEO satellites experienee shorter average eelipse times than EEO 
satellites, thereby simplifying the eleetrieal power requirements of the spaeeeraft. 
Generally, launeh eosts to MEO are slightly more affordable than GEO and allow for 
larger spaeeeraft for the same launeh vehiele. Depending on the speeifie orbit and 
altitude, MEO satellites ean experienee very high proton and eleetron flux within Earth’s 
Van Allen Belts [33]. 

3. Highly Elliptical Orbit (HEO) 

A Highly Elliptieal Orbit (HEO) is an orbit with an eoeentrieity larger than 0.5 
[34, p. 6]. HEO is oommonly used as an intermediate transfer orbit for satellites on their 
way to geosynohronous orbits via Hohmann transfer, but there are some speoialized 
soientifio and eommunieations applioations. During its twelve hour orbit, a HEO satellite 
will dwell for approximately nine hours near apogee before sweeping through perigee in 
the remaining three hours of the orbit. HEO satellites typioally degrade faster than MEO 
or GEO satellites due to the atmospherie drag during their transit near perigee, so some 
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fuel is usually reserved to extend the operational lifetime through regular orbit 
maintenanee maneuvers. Like MEO satellites, HEO satellites ean experienee very high 
proton and eleetron flux within Earth’s Van Allen Belts depending on their speeifie 
orbital profde [33]. 

Eor soientific applications, HEO orbits allow for measuring detailed cross sections 
of interesting phenomena. By transiting across a large range of geocentric distances out 
to twenty Earth radii, these satellites produce volumes of data on the structure and 
composition of Earth’s magnetosphere and the solar-terrestrial interaction. As shown in 
Eigure 21, specially designed 12-hour, “Molniya” orbits are commonly used to provide 
long dwell times over high northern latitudes where the coverage by GEO satellites is 
limited. The same principle could be used over any latitude of specific interest. By 
carefully matching the inclination (± 63.4 degrees) and eccentricity of these orbits, 
engineers establish an orbit that keeps the argument of perigee constant, i.e., the latitude 
at which apogee occurs does not vary. Satellites in Molniya orbits dwell over the 
northern hemisphere during the orbit, allowing for excellent coverage of the northern 
latitudes with fewer satellites in a constellation. As its name suggests, this is a standard 
orbit for Russian communications satellites [35]. 
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4, Geosynchronous Orbit (GEO) 


GEO satellites cirele the Earth at a fixed altitude of 35,786 km and have a period 
of 23.9345 hours, or a sidereal day, matehing the time it takes for the Earth to rotate a full 
360 degrees. The satellite’s velocity at this altitude matches the Earth’s rotational rate, 
thereby keeping the satellite over a specific location on the Earth’s equator and creating a 
satellite ground track in a “figure-eight” pattern on that location. A geosynchronous orbit 
with zero inclination is generally referred to as a geostationary orbit as it will orbit 
directly above a point on the Earth’s equator with no North/South displacement. Many 
satellite communications systems utilize geostationary satellites, which appear stationary 
to a ground observer, to eliminate the need for tracking systems in the ground receivers 
[34]. 

The concept of geosynchronous spacecraft was first proposed by Arthur C. Clark 
in 1945, more than a decade before Sputnik 1, the first satellite, was launched. Clark 
correctly deduced that a system of three satellites in geosynchronous orbit would have 
simultaneous coverage over nearly the Earth’s entire surface. This makes the GEO belt 
an optimum orbit for any application seeking to minimize the number of satellites for 
global coverage; however, there can be significant transmission delays due to the large 
transmission distances involved. 

GEO is highly populated by many communications, weather, and remote sensing 
satellites. Satellites at GEO experience much shorter average eclipse times than satellites 
at lower altitudes and generally avoid the high proton flux existing within Earth’s Outer 
Van Allen Radiation Belts [33]. As discussed in Chapter V, GEO satellites spend most 
of the year in full sunlight and only experience eclipses near the spring and autumn 
equinoxes with each eclipse lasting about 70 minutes, or 5% of the orbital period. GEO 
satellites tend to be very large due to the much greater transmission losses experienced by 
their communications systems at this distance, and their cascading effect on the required 
antenna, electrical, and attitude control system sizes. 

Examples of the orbits discussed above and “typical” satellites within those orbits 
are outlined in Table 5. 
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Orbit 

Low Earth 
Orbit 
(LEO) 

Medium 
Earth Orbit 
(MEO) 

Highly 

Elliptical Orbit 
(HEO) 

Geosynchronous 
Earth Orbit 
(GEO) 

Typical 

Satellite 

Landsat 

7 

GPS 2-27 

Russian 

Communications 

Satellites 

TDRS-7 

Apogee 

altitude 

703 km 

20,314 km 

39,850 km 

35,809 km 

Perigee 

altitude 

701 km 

20,047 km 

500 km 

35,766 km 

Period 

98.8 min 

12 hr 

~12hr 

24 hr 

Inclination 

98.21 deg 

54.2 deg 

63.4 deg 

2.97 deg 

Eccentricity 

0.00010760 

0.00505460 

0.741 

0.000765 


Table 5. Typical Orbit Parameters. After [35, p. 115] 


E. PERTURBATIVE FORCES 

While the original laws discovered by Kepler and Newton accounted for the 
primary acceleration within a satellite’s orbit, there are many other conservative and non¬ 
conservative perturbative forces that affect the motion of an orbiting body. Equations 
(IV. 1) through (IV.6) assumed that the Earth and satellite were spherically symmetric and 
of uniform density and completely ignored the usually minor effects of the other 
perturbative forces. These simplifications are generally valid for most applications where 
a limited number of orbits are involved or only a rough idea of the satellite track is 
required, but for high accuracy applications, such as power beaming, orbital debris 
clearance, or illuminating extraterrestrial bodies, perturbations on orbital motion must be 
taken into account. 

Similarly, it can be an egregious error to assume that perturbative forces are 
always small. Atmospheric drag during reentry and the second body gravitation in the 
restricted three-body problem are prime examples of a perturbative forcing growing to 
match the primary attractive force. Whenever they are this large, they are no longer 
termed “perturbations,” because the fundamental nature of the problem has changed. In 
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this paper, the perturbative aecelerations are discussed because the removal of orbital 
debris, covered in Chapter VI, depends upon increasing the atmospheric drag of the 
debris particles. 

All of the significant accelerations that will affect an orbiting object are discussed 
in the following sections and displayed in Figure 22 on a log-log scale. At the top of the 
plot, the gravitational acceleration due to a uniformly spherical Earth (labeled GM) is 
three orders of magnitude larger than the other accelerations. The J2, J4, and Jg lines 
display the harmonic expansions of the geopotential. The odd zonal harmonics have no 
temporal effect on the Keplerian motion of a satellite, and they are not included here [27, 
p. 92]. J 2 represents the Earth’s equatorial bulge due to its angular momentum. While 
atmospheric drag and lift are virtually non-existent above 1000 kilometers, it can even 
exceed the gravitational potential at altitudes below a few hundred kilometers, allowing 
for winged flight within the atmosphere. As an important note, the gravitational 
attraction due to the Sun, Moon, and other planets produces the only accelerations with 
positive slope, i.e., their forces increase as the distance from the Earth increases. The 
Solar radiation pressure generally drops with distance from the Sun, but on the scales 
shown here, it remains virtually unchanged. Overall, the Tidal, Relativistic, and Albedo 
perturbations remain fairly low and trend together. 
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Figure 22. Central and Perturbative Accelerations. From [27, p. 69] 
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Conservative Perturbations 


a) Attraction by the Earth 

As described by the geopotential expansion shown above in Equation 
(IV.7), the acceleration due to Earth’s gravity is an extremely complex variable to 
determine. Each expansion term adds a level of precision by accounting for increasingly 
subtle differences in the mass distribution of the Earth. The Newtonian central 
acceleration (GM) is inversely proportional to the radius squared, and therefore it is 
represented in the log-log plot by a straight line of slope s=-2. The J 2 acceleration term, 
representing the Earth’s oblateness, is inversely proportional to /, and it is represented in 
the plot as a straight line with slope s=-4. The other high order harmonic terms {Jn) are 
inversely proportional to and they have a slope s=-(n+2). 

b) Attraction by the Sun and Moon 

While a satellite obviously experiences acceleration due to the 
gravitational attraction of the Sun and the Moon, the Earth is also affected by their 
gravity, and therefore, it is the differential attraction felt by the satellite that must be 
accounted for here. Due to its proximity, the Moon exerts the largest, non-terrestrial 
acceleration on the satellite- followed closely by the gravitational acceleration due to the 
Sun’s extremely large mass. In general, the Eunar and Solar differential accelerations 
increase with orbital altitude because the average gravitational attraction to them 
increases. 


c) Attraction by Other Planets 

Planetary perturbations of Earth-orbiting satellites are small and vary with 
the orientation of the planets, with Venus usually providing the largest contribution, and 
they can cause no significant uncertainty in the Keplerian orbit [27, p.71]. At much 
greater ranges than the Moon and lower relative mass than the Sun, the other planets 
exert a differential, gravitational attraction similar to that described above, but it is at 
least 1000 times less than the Solar or Eunar accelerations. 
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d) Tidal Elements 


The gravitational attraction of the Moon, Sun, and other planets also affect 
the Earth’s oceans and crust, causing tidal variations in the liquid and solid mass 
distribution of the Earth. The Earth’s crust can raise and fall with an amplitude of as 
much as one decimeter, affecting the gravitational potential almost ten times more than 
the ocean’s tides [27, p. 71]. The tidal acceleration is significantly greater at lower 
altitudes as the range is lower and the angular differential due to the displaced mass is 
larger. 


e) Relativistic Effects 

While satellites orbit at speeds well below relativistic velocities, Einstein’s 
relativity provides a small correction to the Newtonian equations of motion resulting in 
equally small changes in the Keplerian motion of the spacecraft. Since they are about six 
orders of magnitude less than the J 2 acceleration, these changes are generally negligible 
over a few orbits, but they result in a secular effect on the argument of perigee. Over 
multiple orbits the perigee of an orbit will change more rapidly than classical calculations 
alone would predict. 

2, Non-Conservative Perturbations 


a) Atmospheric Drag 


Eor satellites in EEO, friction with the molecules of residual atmospheric 
gases can provide significant decelerations. These perturbations are extremely difficult to 
model since they require a detailed model of the upper atmosphere, which varies 
according to solar and geomagnetic activity. They also depend heavily on the 
spacecraft’s size, configuration, surface materials, and orientation with respect to the 
direction of travel. The acceleration due to atmospheric drag can be computed using 
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where p is the atmospherie density, v^/ is the velocity of the spacecraft relative to the 
atmosphere, Cd is the dimensionless atmospheric drag coefficient (a function of the 
satellite’s shape), Asc is the spacecraft’s effective projected area, nisc is the spacecraft 
mass, and is a unit vector of the spacecraft velocity relative to the atmosphere. 

b) Radiation Pressure 

The actual radiation pressure experienced by a satellite depends heavily 
upon its size, configuration, and surface materials. A satellite in Earth orbit experiences 
three significant radiation pressures- Solar, Albedo, and Earth’s Infrared. Naturally, there 
is no Solar or Albedo acceleration when the satellite is in eclipse. Solar radiation 
pressure, due to the Sun’s direct illumination of the satellite, is inversely proportional to 
the squared distance to the sun, and it is therefore largely independent of satellite altitude. 
The Albedo radiation pressure, due to the direct reflection of the Sun’s rays off of the 
Earth’s surface, is inversely proportional to the squared orbital radius, but it can vary 
greatly due to the extreme variations in Earth’s surface emissivity. Earth’s Infrared (IR) 
radiation pressure, due to the emission of Earth’s heat in the IR spectrum, varies like the 
Albedo, but it can affect a satellite even in eclipse [27]. 

3, Other Perturbations 

The list above does not constitute an exhaustive listing of the perturbations that 
affect an Earth orbiting satellite, but it represents the largest contributors to variations in 
the predicted Keplerian motion. There are many other physical and theoretical 
perturbations, such as: the thermal re-emission of radiation absorbed by the satellite, 
which usually occurs in an anisotropic way due to the variations in emissivity and surface 
temperature; micrometeorite drag, due to the impact of micrometeorites; the Poynting- 
Robertson effect, due to the fact that even if incident radiation emits isotropically, it only 
does so in the satellite’s reference frame; additional general relativistic perturbations; and 
the theoretical emission of gravitational waves, which is many orders of magnitude less 
significant than the other examples provided [36]. 
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V. POWER BEAMING TO SATELLITES 


While Chapter IV deseribed the orbital motion of a satellite, Chapter V 
investigates the use of a ground-based, high-energy PEL, propagating through the 
atmosphere, to power orbiting satellites. In the design of a satellite, the eleetrieal power 
system typieally reeeives the highest priorities for redundaney, testing, and safety 
margins [31, p. 326]. Of all the systems onboard, it is the first to be plaeed in operation, 
and it remains in serviee for the entire life of the mission. It serves as the literal heart of a 
spaeeeraft, pumping power to every other system and allowing for their operation, and 
when it inevitably fails, the spaeeeraft fails with it. 

The frequent eelipses suffered by every satellite in LEO plaee a very high strain 
on the eleetrieal power system, and they ean drastieally affeet the battery and spaeeeraft 
lifetime. During the design phase, the projeeted number and length of eelipses over a 
satellite’s design life begin a proeess driving the size, weight, eomplexity, and final eost 
of the satellite into an ever inereasing spiral. The concept of beaming power from the 
Earth to a satellite via a high-energy laser has been suggested as an optimal solution to 
minimize the depth of discharge on every orbit, and thereby reduce the battery and solar 
array size, the spacecraft mass, and the overall complexity of the satellite. Unfortunately, 
this concept requires a considerable, and ultimately impractical, ground infrastructure. 
The ground tracks of most satellites vary considerably from pass to pass, and it requires a 
large number of ground stations to ensure reliable power for every nighttime pass. 

Although it may seem desirable, power beaming will not be able to provide the 
constant and economical coverage that is required to meet most of a satellite’s power 
requirements in the foreseeable future. Another futuristic, but more realistic, potential for 
the power beaming concept resides in the ability to improve a satellite’s design lifetime in 
certain specific orbits, such as a sun synchronous orbit. Eor a sun-synchronous orbit, the 
orbital profile is highly regular, and, therefore, it would probably be the least complex 
application of power beaming. Eor this specialized orbit, it might be possible to utilize 
fewer ground stations, thereby reducing the overall investment in ground infrastructure. 
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By reducing the battery’s depth of discharge in every orbit, this sun-synchronous, 
power-beaming system could allow for a smaller and lighter electrical power system. 
The cost savings for extending the lifetime of an individual spacecraft can be significant, 
but additional cost benefits could be realized by designing smaller and lighter spacecraft 
and using a less capable and cheaper launch vehicle. Current launch costs to LEO are 
approximately $8,400/kg [44, pg. 178]. If consistent power beaming allowed for a 
satellite design with 10% less mass, then a 450 kg satellite might save about $833,000 in 
launch costs. Despite these considerable savings for an individual spacecraft, the large 
ground infrastructure required for effective power beaming can only become economical 
when applied to satellites on a grand scale. Even for a typical EEO satellite constellation 
of ten-to-twenty spacecraft, the cost of infrastructure required to support a worldwide 
capability for power beaming is prohibitive. So, it may only be cost effective to support a 
full ground based infrastructure as a national or corporate enterprise (by servicing more 
than one constellation of satellites). 

Currently, the U.S. Navy has funded the development of a 100 kW class EEE with 
intentions to pursue the creation and deployment of a MW class laser for shipboard use 
[37]. Presumably, the lasers onboard these ships would remain deployed and operational, 
and for the majority of the time, their lasers would not be occupied by combat operations. 
Although a ship in a maritime environment would operate under less than ideal 
conditions for long-range, atmospheric propagation, a deployed fleet of naval vessels 
with MW-class EEEs might be the preferred terrestrial infrastructure for powering critical 
government, military, or communications satellites in specialized orbits or during 
emergent situations. Eor military or national security applications, payloads must 
typically restrict their operations during eclipse due to their limited battery power. The 
ability to power critical missions during eclipse might have a considerable effect on 
battle-space awareness or control. In special cases, the additional power during eclipse 
could allow for critical communications during wartime or additional imaging by remote 
sensing satellites. In emergency cases, power beaming may even assist in the recovery of 
satellites that have lost attitude control, failed to fully deploy their solar panels, or are 
otherwise unable to maintain their solar panels in a sun-soaking orientation. Earth-based 
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lasers could provide additional irradiance to the spacecraft, helping to provide useful 
power and heat while engineers attempt to recover the craft. 


A, ECLIPSE DISCUSSION 

The fraction of time that a satellite spends in sunlight and eclipse is important to 
the design of the electrical power and thermal systems. Depending on the time of year, 
orbit, and spacecraft age, there can be significant differences in a satellite’s eclipse 
profde and power budget, so a spacecraft is typically designed for maximum eclipse and 
minimum sunlight at the end-of-life (EOL). Thus, every satellite is over designed for all 
other conditions within its mission lifetime. Similarly, the orbital and eclipse periods 
determine the total number of battery cycles and depth of discharge per orbit, which are 
directly related to the battery lifetime. 

The typical eclipse profile of a low earth orbiting satellite in a circular orbit at 400 
km altitude and zero degrees inclination is shown in Figure 23. This plot displays the 
angle from the orbital plane to the sun (Ps), or “beta angle,” overlaid with the angular 
radius of the Earth from the altitude of the satellite’s orbit. For an orbit at zero degrees 
inclination, the beta angle will merely vary with the seasonal variation of the Sun’s 
declination (± 23.5 degrees), and at 400 km, the Earth’s angular radius is approximately 
70 degrees. When the beta angle is inside the Earth’s angular radius (plotted in red), the 
satellite will pass through the Earth’s shadow during each orbit, and the absolute 
difference between the values is proportional to the length of each eclipse. When the beta 
angle is outside of the Earth’s angular radius, the satellite will not enter Earth’s shadow, 
remaining in full sunlight for every orbit. As shown in Figure 23, this EEO satellite will 
pass through eclipse on every orbit throughout the year although at some points of the 
year the time in eclipse will be less than others. 
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In Figure 24, the typical eclipse profile of a LEO satellite in a circular orbit at 400 
km altitude and 65 degrees inclination is plotted. In this plot, the precession of perigee 
for the orbit combines with the seasonal variation in the Sun’s declination to produce 
about 30 days without eclipse for the satellite, and for some of the year, the satellite 
experiences much shorter eclipses than those faced by the LEO satellite with zero 
inclination. 

In contrast to the eclipses experienced by LEO spacecraft, the graph of a typical 
geosynchronous spacecraft shown in Eigure 25 shows that it will only suffer eclipses 
during about three months of the year, and these eclipses are very short compared to the 
overall orbital period. Since most GEO satellites orbit near zero inclination, the beta 
angle varies with the seasonal variation in the Sun’s declination, but at GEO altitudes, the 
angular radius of the Earth is only about 8 degrees. 
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Beta vs Day of the Year 



Figure 24. LEO Eclipse Profile 65 Degree Inclination. After [38] 


Beta vs Day of the Year 



Eigure 25. 


GEO Eclipse Profile. After [38] 
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Figure 26 shows the proportional time in eclipse for each orbit discussed above. 
In a LEO spacecraft, eclipse can be as much as 36 minutes out of a 90-minute orbit 
(eclipse for 39% of each orbit), and, depending on the mission, power budget, and battery 
type, the spacecraft may not recover a complete charge prior to the next eclipse. For a 
GEO spacecraft, the eclipse can last as long as 70 minutes out of a 1436-minute orbit 
(eclipse for ~5% of each orbit), and the satellite can easily recover a complete battery 
charge prior to the next eclipse. In both cases, power beaming could limit the cyclical 
stress on the electrical power system, but each orbit type has individual advantages. With 
eclipse times almost twice as long as a EEO satellite, a GEO satellite with a similar 
mission and loads would need a much larger battery. 
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Figure 26. Proportional Time in Eclipse for EEO and GEO Satellites 
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B, POWER BUDGET 

When designing a spaeeeraft and planning for on-orbit operations, all of the systems 
and eomponents must be designed to the most limiting ease and still be able to meet mission 
requirements. Sinee every spaeeeraft is adversely aflfeeted by radiation, eebpses, and other 
faetors while in orbit, this means that the typieal spaeeeraft is launehed with enough exeess 
eapaeity at beginning-of-life (BOL) to ensure that it ean meet all mission requirements at 
end-of-life (EOL), when the solar eell and battery effieieneies have degraded. 

Shown in Figure 27 is an approximate power budget, or skyline ehart, for an imaging 
satellite, displaying the expeeted systems in operation and power required for eaeh system 
over one orbit. The peaks in Figure 27 ean be direetly related to battery eapaeity and solar 
eell sizing as these eomponents must supply the power to meet all of the requirements plaeed 
upon them during eelipse or sunlight. The average spaeeeraft spends its time in the sun 
operating the payload, oommunieating with ground stations, and power soaking the solar 
arrays to eharge the battery in preparation for the antieipated battery diseharge during eelipse. 
During this power soaking, or sun soaking, of the solar arrays, the satellite maintains an 
orientation that maximizes the solar illumination of the solar arrays. While ehanging 
orientation for oommunieations, imaging, or sun soaking, the satellite’s attitude eontrol 
system utilizes power to point the spaeeeraft. During eelipse, the satellite’s thermal eontrol 
system operates extra heaters to stabilize the temperature of the satellite. 

Obviously, the design of the eleetrieal power system and the availability of solar 
power ean eonsiderably limit the operation of a satellite. Remote sensing satellites typieally 
have a restrieted window for operating the payload during the daytime, and nighttime 
operations are tightly eontrolled due to the limited battery eapaeity. By eontinuing the 
battery eharge after a satellite enters eelipse, a laser eould potentially improve eurrent satellite 
operations. The primary differenee between the imaging satellite depieted in this figure and a 
eommunieations satellite is that a eommunieations satellite often has a nearly eonstant power 
profile due to the requirement for eontinuous oommunieations ability, i.e., eonstant operation 
of the payload and attitude eontrol systems. 
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Figure 27. Generic Power Profile of an Imaging Satellite. After [38] 


As mentioned above, the total battery capacity becomes a dominant parameter 
early in the design phase of a spacecraft, as it must accommodate the power requirements 
for every component for the entire spacecraft lifetime. In Figure 28, the battery capacity 
for a typical spacecraft is plotted versus time to display the frequent charge and discharge 
periods that a satellite battery will undergo. A satellite is not always able to fully 
recharge its battery to maximum capacity during each orbit. Mission requirements 
sometimes dictate more frequent payload operation and less time for sun soaking, 
resulting in an incomplete charge while in sunlight. Likewise, in a process called 
“battery conditioning,” engineers will increase the depth of discharge over a few orbits to 
preserve the overall battery capacity. Historical failure data suggests that the depth of 
discharge per orbit and the total number of battery cycles are quite important to the 
overall lifetime of a spacecraft as they are directly linked to the stresses on the electrical 
power system. 
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Figure 28. Typical Battery Capacity Graph. From [38] 

C. PHOTOVOLTAIC CELLS 

By far the most common method of providing power to a satellite is through the 
use of solar panels, which utilize large arrays of specialized photovoltaic cells to gather 
some of the energy provided by the sun’s illumination. Photovoltaic cells utilize the 
photoelectric effect to transform incident electromagnetic radiation into electrical current. 
Solar cells are photovoltaic cells specially designed to utilize incident solar energy for 
power, and increasingly, multi-junction solar cells are becoming common in satellite 
applications. Multi-junction cells attempt to capture more of the energy within the broad 
solar spectrum by mixing materials with differing band gaps. Regardless of their design, 
most of the radiation that falls upon a solar cell will always be wasted due to the breadth 
of the solar spectrum and the limited band gaps of the semiconductor materials. Since 
their introduction early in the space age, solar cell efficiencies have remained below 20% 
for most applications with typical values around 15% [39, p. 277]. 

When designing a photovoltaic cell for laser power beaming, it would be possible 

to better match the band gap of the cell with the wavelength of the incident laser beam, 
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thereby redueing exeess energy delivered as heat and inereasing the overall effieieney of 
the eell. Photovoltaie eells operating at 800 nm ean be subjeeted to more than five times 
the natural solar radiation intensity (-1366 W/m ) without requiring aetive eooling, and 
more than 60% effieieney has been demonstrated using Gallium Arsenide (GaAs) eells at 
this wavelength [39, p. 277]. Generally, the eonversion effieieney of a speeially designed 
photovoltaie eell under monoehromatie illumination ean be two-to-three times greater 
than a typieal solar eell under solar radiation, but even eonventional solar eells ean 
benefit from monoehromatie illumination. Figure 29 shows a graphie of observed solar 
eell effieieneies when subjeeted to monoehromatie illumination. 



Wivl^ogth, \iMn 

Figure 29. Solar Cell Effieieneies Under Monoehromatie Illumination. From [40, p. 

235] 
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Power beaming has the potential to greatly enhanee the power generated by a 
satellite without inereasing the size of the solar arrays. Using the Hubble Spaee teleseope 
as a typieal imaging satellite, a rough eomparison of the power generated by solar and 
laser illumination ean be performed. The average solar radiation intensity in Earth orbit 
is approximately 1366 W/m , and the Hubble solar panels oeeupy a total area of about 
38m [41]. Thus, the total power falling on these panels is 52 kW, but even with a 
generous effieieney of 20%, the solar panels would only provide 10.4 kW of eleetrieity. 

By applying the high intensities allowed through monoehromatie illumination, a 
laser eould eonservatively projeet more than three times the solar intensity to the solar 
panels, about 4200 W/m , without requiring additional radiators. On the same area, this 
illumination results in roughly 160 kW falling upon the solar panels, and with an 
estimated monoehromatie effieieney of 50%, the solar panels would generate about 80 
kW of eleetrieity. Under monoehromatie illumination, the eombination of higher eell 
effieieneies and higher sustainable radiative intensities allows the satellite to generate 
almost 8 times more power than it eould under solar radiation. This additional power 
eould be erueial for speeialized applieations sueh as providing a power “boost” to eounter 
jamming of GPS or oommunieations satellites or rapid eharging of a satellite’s battery 
during eelipse. 

Solar eells experienee eonsiderable damage due to the eleetron and proton flux in 
Earth orbit and the thermal stresses of frequent heating and eooling. This damage ean 
seriously degrade the EOE effieieney and power output of the solar panels, neeessitating 
larger design margins and solar panels at launeh. Some seientists have suggested 
illuminating solar panels with a high power laser to anneal the eells and repair some of 
this damage remotely. In this instanee, the targeting aeeuraey of the laser mount and the 
thermal effeets on the rest of the spaeeeraft body are mueh more prominent eoneerns than 
for power beaming. Otherwise, this applieation is similar to power beaming with mueh 
less frequent aeeess to the spaeeeraft, i.e., it would be performed infrequently. Thus, only 
one or two ground sites would need a laser eapable of the higher power and better 
aeeuraey that is required to support the “repair” illumination. 
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D, 


SOLAR RADIATION SPECTRUM 


Although the sun outputs a great amount of energy, it does so aeross the entire 
eleetromagnetie speetrum. In eontrast, the PEL ean transmit its energy within a 
bandwidth of about 0.1% of the wavelength, or about 10' pm for an FEE operating at 1 
pm [15]. The radiation output of the sun can be easily approximated using the Stefan- 
Boltzman law by assuming that the sun radiates as a black-body radiator at a temperature 
of about 5800 K. This temperature approximates the temperatures of the outer visible 
layers of the sun’s corona, those that directly interact with the solar system through 
radiative means. The Stefan-Bolztman law is 


Rsun = , (V.l) 

2 8 
where Rsun is the power radiated in W/m , a is the Stefan-Boltzman constant (5.67 x 10' 

W/m^K"^), E is the emissivity (approximately one for a blackbody), and T is the 

temperature in K. Using Equation (V.l) and the assumptions above, the total output 

power of the sun is approximately 64.2 MW/m at the solar surface. For the applications 

discussed here, it may be more useful to know the radiative flux as a function of 

wavelength since the laser only has to compete in limited bandwidths. [33] 

Planck derived an analytical expression for the radiative flux as a function of 
wavelength for a blackbody. 


R 


SunX 



(V.2) 


2 

where Rsuni is the blackbody thermal radiation in W/m , X is the wavelength in m, C/ is 
3.74 X 10'^^ kW-m^, C 2 is 1.439 x 10'^ m-K, and T is the temperature in K [42]. This 
relation can be used to generate the radiation curve of a blackbody, shown in Figure 30, 
but it can also be used to determine the sun’s power over a limited bandwidth. In most 
applications, the laser only has to compete with the sun within its bandwidth or near the 
band-gap energy of a photovoltaic cell. 
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Figure 30 above compares the actual irradiance of the sun in Earth orbit, referred 
to as Air Mass Zero (AMO), over a range of wavelengths to the power profile of the 
theoretical blackbody radiator. It also shows the ground level irradiance in red, 
highlighting the absorption caused by Earth’s atmosphere on the AMO radiation. As 
shown above, the actual irradiance is very similar to the blackbody approximation at the 
wavelengths of interest for high energy applications and atmospheric propagation. 
Calculations in Chapter VII will utilize the blackbody approximation to determine the 
sun’s output within the laser bandwidth. While satellites in low earth orbit do operate 
inside the upper vestiges of the atmosphere, atmospheric effects on incident radiation are 
generally negligible, and the AMO solar spectrum can be assumed for most spacecraft. 
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E, 


ATMOSPHERIC PROPAGATION 


When attempting the beaming of power from a terrestrial site to an orbiting 
spaeeeraft, the absorptive and seattering properties of atmospherie water moleeules begin 
to oritieally affeet laser propagation. As shown above in Figure 16 and Figure 30, the 
total transmittanee of Earth’s atmosphere near 1 pm ranges from 50 to 70% of the total 
transmitted power. Thus for the simulations and ealeulations in Chapters V, VI, and VII, 
we must assume that whatever the final output of the laser at the direetor mirror, only a 
eonservative 50% will pass through the atmosphere. Still, when we eompare the 
antieipated 1 MW output power of a ship-home FEE to the typieal tens of kilowatts 
required for even the largest spaeeeraft, this loss will not prevent sueeessful power 
beaming. 

Eor GEO applieations, the laser beam will remain relatively motionless in the 
atmosphere and eould experienee thermal blooming. In the EEO ease, the angular rate of 
the laser beam should avoid thermal blooming high in the atmosphere, but thermal 
blooming eould develop in the dense and relatively stagnant sir near the direetor. Some 
method of indueing eonveetion near the beam may allow for a reduetion in thermal 
blooming, sueh as the use of large turbofans to move air through the laser eolumn. This 
artifieially indueed “wind” would tend to bend the laser beam whieh eould possibly be 
eorreeted by realigning the beam direetor. 

F, BEAM SPOT SIZE AND INTENSITY 

Satellites speeially designed for power beaming applieations might utilize a single 
eireular photovoltaie panel to assist in laser targeting and reduee illumination of the 
satellite body, but most eurrently orbiting satellites utilize dual, sun-pointing solar arrays 
extending symmetrieally from opposite sides of the satellite. Therefore, any power 
beaming to a generie satellite will have to individually target eaeh solar panel or generate 
a spot size large enough to illuminate the entire spaeeeraft body and its solar panels. For 
a generie, imaging satellite like the Hubble Spaee Teleseope, the total diameter of the 
beam required to eneompass both solar panels is about 13 m [41]. 
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Using Equations (III. 1-6), the on-orbit beam radius and intensity of a generie, 
MW elass, ship-borne laser ean be estimated without turbulenee. The beam radius (w) is 
given by 


w{z) = 



, and 


(V.3) 


the intensity (7) is given by 



TTW^ ■ 


(V.4) 


The orbital altitude of a satellite gives the minimum propagation range for a laser 
engagement, whieh is assumed to be 400 km for the simulations in this ehapter. Figure 
67 in Appendix B shows the maximum range of engagement as 1844 km. The beam 
direetor, whieh determines the beam waist radius {wo), is assumed to be approximately 1 
m. This results in a Rayleigh range {zr} of 3.14 x 10^ m. For the maximum FEO range 
of 1844 km, the ealeulated beam radius at the target is 1.16 m, and for the minimum FEO 
range of 400 km, the ealeulated beam radius at the target is 1.01 m. Atmospherie 
transmission is assumed to be a eonservative 50 % of the transmitted power, and the 1 
MW of transmitted energy results in about 500 kW of power delivered to the satellite. 
By diffraetion alone, the respeetive beam intensities for the maximum and minimum 
range engagements are about 118 kW/m and 156 kW/m . A eonservative estimation of 
atmospherie turbulenee and thermal blooming might halve these intensities, but this still 
results in about 100 kW/m ineident on the satellite body and solar panels that are 
designed to reeeive solar irradianees of ~1.4 kW/m [15]. 

These intensity values exeeed the power requirements of most orbiting satellites 
by a few orders of magnitude, and they eould result in damage to the satellite. Therefore, 
for power beaming applieations, the laser power may need to be lowered or the beam 
direetor may need to be redesigned to deerease the beam intensity at the satellite. 
Fortunately, the FEE is sealeable over a wide power range, and they eould be designed to 
“dial-in” powers between 1 kW and 1 MW [15]. Fikewise, the beam direetor is already 
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expected to be curved and utilize an adaptive optics system to improve atmospheric 
propagation effects. On days when turbulence significantly reduces intensity, the MW 
level PEL can still be effective. 

For satellites in GEO, the propagation range is about 35,678 km. By diffraction 
alone, the calculated beam radius at the target satellite would be about 11.4 m, and the 
beam intensity is about 1.22 kW/m^. This is slightly less than the solar intensity, but it 
may be enough to power vital systems throughout eclipse for most satellites. However, 
in this case, thermal blooming effects occur much closer to the laser source than to the 
target, and over these long propagation ranges, they can have a much greater effect on the 
beam intensity at the target. The combined effects of thermal blooming and turbulence 
might make power beaming to GEO impractical for beam directors of this size. If the 
beam director mirror’s radius is doubled, then the beam radius at the satellite is halved, 
and the beam intensity is about 4 times larger than the solar intensity. 

In the Hubble example discussed above, a deliberate spreading of the beam to an 
on-target diameter of 13 m would reduce the overall intensity to about 3.77 kW/m , and 
further atmospheric effects may lead to on-orbit intensities of about 2 kW/m . This is 
only slightly more than the current incident solar irradiation and should be well within the 
satellite’s design tolerances. However, this approach utilizes basic assumptions about 
atmospheric transmission and a simplistic illumination pattern of the target satellite. Any 
specific application of power beaming to an orbiting satellite would have to accurately 
measure and anticipate the atmospheric propagation and model the beam intensity 
fluctuations across the satellite body to prevent unintentional damage. A dedicated site 
using a MW-class FEE might use a larger mirror to accommodate smaller spot sizes, and 
therefore, larger intensities at distant satellites. 

G. SATELLITE THERMAL EFFECTS 

Currently, the thermal control system of every satellite is designed to withstand full 
solar illumination through much of their orbits and possibly even long periods of sun-soaking 
in a fixed orientation. Thus, the coatings and materials of every satellite surface are carefully 
analyzed and controlled. Assuming this full solar illumination, the major constraints on the 
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solar panels are that they have a large enough area to support the satellite loads, battery 
eharging, and some method of aligning themselves to the solar radiation. For a satellite 
designed for power beaming, a single cireular solar panel would be the optimal target for the 
laser illumination, while still remaining available for the eolleetion of solar illumination. 
There will be changes to the thermal control system, coatings, and materials of future 
satellites designed for monochromatic, laser illumination. 

The radiative heat sources contributing to a satellite’s final temperature are 
graphically depicted in Figure 31 and included in Equation (V.5). The total heat input to 
a typical satellite is comprised of terms representing the direct solar radiation, albedo (a 
reflection of solar energy directly off the Earth), Earth’s infrared (IR) radiation, 
equipment heat dissipation, and solar cell efficiency (since the solar cells turn some 
incident energy into electricity not heat). The heat input equation is 

Input SolarRad Albedo EarthlR Equip SolarCell 

Qi„=ot^s^p + PaaIsAF,+eIjj,AF^ + Q,^-rH , (v.s) 

where Qm is the total heat input to the satellite, a is the absorptivity, 4 is the solar 
irradiance, Ap is the projected area, pa is the albedo factor, A is the surface area. Fa is the 
Albedo view factor, s is the emissivity of the surface, Im is the Earth IR irradiance, Fg is 
the Earth view factor, Qgq is the heat input due to operating equipment, rj is the solar cell 
efficiency, and fp is the cell packing factor. In this case, albedo refers to the direct 
reflection of solar energy off of the Earth’s surface, which is differentiated from the 
general thermal IR radiation of the Earth’s IR term [38]. 
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Figure 31. Satellite Thermal Interactions with the Space Environment 


Power beaming during daylight could be useful in propulsion and other 
specialized applications. If power beaming is conducted while the satellite is in sunlight, 
then the solar cells would be pointed towards the Earth and operate with a different 
efficiency under monochromatic illumination. In this case, there is an additional term 
representing the heat deposited by the laser on the spacecraft and the final term is 
replaced by the laser energy transformed into electricity by the photovoltaic cells. The 
heat input equation for daylight power beaming is then 




SolarRad LaserRad Albedo 

ol I „ -\- ot IT Q cx I F 

s p P r a s a 


EarthIR Equip PVCell 

+el jr AF^ + Q^q -Vpv^L^pfp 


(V. 6) 


where II is the laser irradiance at the target and rjpv is the photovoltaic cell efficiency. 
Power beaming while the satellite is in the sun is a much more complex thermal problem 
as the satellite will be subjected to illumination due to the sun, laser, albedo, and earth’s 
IR from various angles. As shown in the equations above, thermal control of the 
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spacecraft depends on the particular absorptivities and thermal conductive properties of 
the satellite materials, but generally, power beaming to a satellite in sunlight will add 
significantly to the heat load on the satellite. Therefore, power beaming during sunlight 
may only be worthwhile for applications that require more power than solar illumination 
alone can provide, such as for propulsion. Likewise, it is probably feasible only for 
specially designed satellites with appropriate photovoltaic panels, surface coatings, and 
possibly extra radiators to dissipate the additional heat deposited by the laser. 

The thermal profile of a satellite receiving laser illumination during eclipse would 
have terms similar to those in Equation (V.5). This equation replaces the solar radiation 
term with a lasing radiation term, removes the albedo term, and adds a heater term to 
account for the heat required to maintain safe operating temperatures for equipment. The 
heat input equation for power beaming during eclipse is 

Input LaserRad EarthIR PVCell Heater 

Qin = + ^IlR +Qeq-Vpv^L ^pfp + Qhtr ’ 

where Qhtr is the heat input by the spacecraft heaters. In this case, the high efficiency 
photovoltaic cells transform most of the incident radiation that falls on them into 
electricity, but some of the remaining energy is deposited as heat in the solar panels. 
Likewise, the beam spot size on-orbit may encompass much of the satellite body, and 
some of this energy will also be deposited as heat. Certainly, satellite surfaces could be 
specially designed for higher or lower absorption near the laser wavelengths, but this 
would require a detailed analysis of the satellite’s thermal profile. Reflecting a specific 
known wavelength from the laser can remove much of the heat absorbed on parts of the 
satellite other than the photovoltaic panels. 

H. PROPULSION 

The current average estimate of launching a satellite to LEO is about $8,400/kg, 

and the cost for the transfer to GEO using a typical, Hohmann transfer is over $21,000/kg 

[44, p. 178]. A Hohmann transfer is a method for changing orbital altitude that utilizes 

two bums on opposing sides of the orbit. Astoundingly, the fuel for a traditional orbital 

transfer can represent over 50% of the total launch mass of the satellite, which is a 
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significant use of money and launeh assets. In an attempt to lower the traditional eost 
and risks of launehing spaeeeraft into orbit, NASA has reeently funded a spaee elevator 
eontest to reward the engineering development and practieal demonstration of the 
teehnologies required to ereate a robotic climbing vehicle that would ascend an extremely 
long eable into spaee while powered by a high energy laser. Some NASA seientists think 
that it may also be possible to eonduet a satellite’s transfer from LEO to GEO using a 
high energy laser to beam power to a “space tug” with an ion thruster engine [45]. This 
space tug would remain in orbit, reduce launeh costs, reduce launch weight, and deliver 
vehicles to GEO with much more payload mass or station-keeping fuel. 

Eigure 32 displays the eross seetions of the EEO to GEO orbital transfer 
maneuver using a traditional, (a) Hohmann transfer or (b) the laser-powered, spiral 
transfer. Unlike the Hohmann transfer, whieh is traditionally eondueted using a liquid 
fueled engine for two large burns, the Spiral orbital transfer maneuver using an ion 
thruster provides exceptionally high efficiencies because the specifie impulse (Lp) is 
large. The fip is proportional to the change in momentum per unit amount of propellant 
used. The primary advantage of using laser versus solar energy for powering the orbital 
transfer is in the mueh larger irradiances and efficiencies that it can provide in lower 
orbits and the nearly eonstant power available near GEO. This would allow for a more 
rapid trip from EEO to GEO, thereby limiting the overall exposure time of the satellite to 
the high proton and eleetron flux within Earth’s Van Allen Belts and putting the asset 
into serviee more rapidly. Easer powered orbital transfers to GEO eould last less than a 
month, whieh is a mueh better option than the several month voyage and high radiation 
levels of an equivalent solar powered system [46, p. 24]. 
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Figure 32. Hohmann and Spiral Orbital Transfers. From [46, p. 26] 


Unlike eonventional missions whieh operate from existing launeh pads and 
benefit from an established ground infrastrueture, laser powered propulsion would 
require a signifieant initial investment in the ground infrastrueture. However, as the 
number of laser powered missions grows, the eost of the infrastrueture is spread among 
more and more launehes and the eost per mission deereases eonsiderably. As shown in 
Figure 33, a laser powered vehiele offers a mueh higher effeetive fp than eonventional 
vehieles as well as signiheant savings on the per kg eost of transportation to GEO. Like 
other spaeeeraft driven by eleetrie propulsion, a laser powered vehiele would have a high 
effeetive fp, but the laser powered vehiele eould presumably use a mueh smaller and less 
massive power plant than solar-powered, eleetrieahy-driven spaeeeraft. Therefore, the 
less-massive, laser-powered spaeeeraft ean have a mueh larger Av and transit through the 
Van Allen Belts more rapidly than other spaeeeraft using eleetrie propulsion. 
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• FOUR GROUND SHTES 

• PAYLOADS: 20/YR (O 2500 kg) 

• EARTHIEO LAUNCH COST s $11 K/kg 


• LASER POWER VEHICLE^ Isp a 2500 

• CONVENTIONAL VEHICLE^ Isp > 480 
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Figure 33. GEO Transportation Costs. From [46, p. 24] 


In one major advantage over conventional launches, the space tug remains useful 
even after its satellite is delivered to GEO. After undocking, the space tug could be used 
to de-orbit other satellites at EOF or to collect and de-orbit large pieces of orbital debris. 
This process could be repeated until the electric drive runs out of propellant. As 
discussed in Chapter VI, debris in high Earth orbits can remain in orbit for millennia, and 
it presents a hazard to all future space missions in crossing orbits. At the present time, 
there are no active methods available for the removal of space debris, only recent 
procedural changes to minimize their future creation, so this use alone may justify the 
expense of designing and building a large, terrestrial, power-beaming infrastructure. 

I. CONCEPT OF OPERATIONS 

There are many claims within the scientific and engineering communities about 
the potential applications of power beaming technology. Most of these applications, such 
as eliminating satellite batteries and powering them through every eclipse, require 
systems of perfect reliability and immense complexity, and they will remain beyond our 
capabilities for many years to come. However, there are some applications that would be 
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well within our technological capabilities assuming the deployment of a MW class laser, 
and many that could even function with modern lasers in the hundreds of kW range. The 
primary challenges facing the implementation of these technologies will be the cost, 
complexity, and safety considerations of building and operating the ground infrastructure 
for power beaming. In addition to the difficulties of finding a ground site that offers clear 
viewing conditions, good atmospheric transmission, and a nearby power source capable 
of supporting a CW, high-power laser, there are significant safety considerations in 
clearing the local area and airspace prior to firing a high energy laser. 

Power beaming, if successfully implemented, offers the potential to improve the 
entire power budget of the spacecraft during the orbital passes when it receives power 
during eclipse. Due to the restrictions on the satellite power budget, daytime operations 
are sometimes limited to allow for enough time to “sun soak” the solar panels and 
achieve enough battery charge to last through eclipse. By providing illumination during 
eclipse, the additional power not only aids the nighttime power budget, but it can also 
allow for a more aggressive usage of the payload during the daylight pass just prior to the 
over flight of the ship-borne laser. 

The United States Navy has recently funded the development of an Innovative 
Navy Prototype for the MW class PEL, and these systems may be deployed to naval 
surface combatants within the next twenty or thirty years. While these lasers would 
operate at sea level in a maritime environment that is not optimal for beam propagation, 
they would presumably offer the benefits of reliable power, rapid and accurate beam 
directors, a proven sensor capable of tracking objects in orbit, and mobility. Assuming 
the successful design, testing, and deployment of these lasers, this thesis models a fleet of 
up to four ship-borne FELs for power beaming simulations in the Satellite Tool Kit 
(STK) program. The laser will be assumed to operate at 1 pm with the ability to access 
any satellite that is at least 10 degrees above the horizon. The lasers will be located at 
sites that will maximize satellite access in an attempt to “bound the problem” and deduce 
which applications are most useful or achievable. 
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J. 


SIMULATIONS 


Three simulations are illustrated here to determine the useful boundaries and 
teehnologieal eapabilities of near-term power beaming applieations. For all simulations, 
the “ground stations” are assumed to be naval vessels equipped with a high energy laser, 
integrated sensors, and aeeurate traeking systems, and their loeations are ehosen for 
general isolation from inhabited areas and reliable aecess to the satellite within eelipse. 
The simulations examine the number and length of nighttime aeeesses to the satellite to 
determine whether signifieant power levels ean be delivered to a satellite for speeialized 
operations during eelipse. Laser illumination on the satellite while it is in the daylight 
may require satellites with speeial eoatings or additional radiators to safely distribute the 
additional heat load, so daytime illumination is not eonsidered here. These simulations 
are intended to be applieable to any eurrently deployed satellite upon the development 
and deployment of the appropriate lasers to support power beaming. 

Other than by loeating the lasers in the open oeean, these simulations do not 
aeeount for the obvious eomplieations of safely firing a high power laser above the 
horizon. Direet illumination by a laser beam, or even the refleetion or seattering of sueh 
light, eould be highly dangerous for satellites, aireraft, ships, and personnel, and some 
method of safely “elearing the range” must be established prior to the implementation of 
these teehnologies. 

1. Simulation One: The Sun Synchronous Satellite 

Simulation One represents the simplest useful ease for power beaming, a satellite 
in a sun synehronous orbit that repeats its ground traek daily at the same loeal times. For 
this orbit, power beaming eould be attempted with a single ground station plaeed in an 
advantageous geographie loeation beneath the ground traek where the satellite is in 
eelipse. This ground station should have reliable and repeating daily aeeess to the 
satellite while it is in eelipse, but would not have aeeess to the satellite on every eelipse. 
At most, by seleeting a ground station loeated between two adjaeent ground tracks, the 
laser may be able to power the satellite on two successive passes through eclipse in a 

critical mission area, i.e., just prior to passing over Asia from the Western Pacific Ocean. 
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Figure 34 displays the two-dimensional satellite ground track for Simulation One 
in orange and the access between the ship and the satellite in white. For this simulation, a 
sun synchronous satellite at 400km altitude, labeled Imagerl, with a repeating ground 
track was created using the STK Orbit Wizard function. The ship, labeled DDL, began 
the simulation beneath the satellite ground track in the South China Sea, and during the 
simulation, the ship transits a repetitive track beneath the satellite’s nightly over-flight. 
Specific orbital parameters, ship route, simulation settings, and detailed simulation results 
are provided in Appendix B. 



Conceptually, the ship is assumed to have access to accurate satellite ephemeris 
data from an external source in order to acquire and track the satellite in the limited time 
that the satellite is above the horizon. Thus, the access between ship and satellite was 
required to be at least ten degrees above the horizon in order to allow for the ship’s 
sensors to gain track and for range safety to avoid possible laser scattering effects. 
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Similarly, the plotted accesses in this simulation pass over densely populated portions of 
Indonesia and mainland China. The operation of this laser might be prohibited except in 
times of war or broad international cooperation due to the hazards it may present to 
foreign property or populace. 

Figure 35 displays the computed satellite access times between the ship and 
satellite over a one week timeline. As expected for a sun synchronous satellite and a 
relatively stationary “ground” station, the accesses occur at virtually the same local time 
every day, while the satellite is in eclipse, and they last for just over 8 minutes each. 
Thus, for a satellite that experiences just about 36 minutes of eclipse each orbit, the ship 
can power the satellite for more about 22% of the eclipse period during one orbital pass 
per day, out of the nearly 16 orbital passes per day. As discussed above in Section C, the 
satellite may be able to generate power levels as much as 8 times greater than its daytime 
levels, allowing for continued imaging, communications, or other specialized power 
boosted modes of operation. 
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Figure 35. STK Ship to Satellite Access Report for Simulation One 
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Simulation Two: Double Access to a Sun Synchronous Satellite 


Simulation Two re-creates the same sun synchronous satellite used in Simulation 
One, but the laser-equipped ship follows a course parallel to the satellite ground tracks to 
obtain access on two successive nighttime orbits. While Simulation One determined an 
approximate maximum access for a satellite under ideal conditions, Simulation Two 
balances the access times on successive orbits to increase the number of nighttime 
accesses. In this case, the individual access times are slightly less, but this maximizes the 
ability to fully utilize a satellite during eclipse in one specific geographic region, such as 
successive orbits over the Korean peninsula and the Chinese mainland. 

Figure 36 displays the two-dimensional ground track of the sun synchronous 
satellite, Imager2, labeled in blue, and the track of the laser equipped ship, DDL2 in 
yellow. Much like in the first simulation, the actual path of the laser on each access 
would pass over the inhabited islands of Indonesia and might violate foreign airspace. 
This presents a significant challenge to the design and implementation of power beaming 
technology since applications on this scale will certainly require worldwide or at least 
regional cooperation among our allies. In this case, the same results might be obtained 
with a ground station in Australia where the central portion of the country is largely 
uninhabited and airspace control might be achievable. 
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Figure 36. Two-dimensional Ship Position, Satellite Ground Track, and Accesses for 

Simulation Two 


Figure 37 displays the ship to satellite access times over a one week period. As 
expected, the ship gains access to the satellite twice per day on successive passes with 
each access lasting just over five minutes. Thus, the satellite only receives about 62.5% 
of the illumination time provided in each pass by the ideal case in Simulation One, but it 
still remains illuminated for about 14% of its time in eclipse. Depending on the specific 
mission, this additional five minutes at the higher powers provided by a monochromatic 
laser could allow for extra payload operation during eclipse or make up for some extra 
usage during its time in the sunlight. 
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Figure 37. STK Ship to Satellite Aeeess Report for Simulation Two 


3, Simulation Three: The Mid-Inclination Satellite 

Simulation Three simulates a satellite in an orbit with moderate inclination and a 
variable ground track. This scenario requires more laser “ground” stations and offers less 
reliable access to the satellite, but these orbits are useful for satellites that are dedicated to 
global coverage vice repeated surveillance of a particular geographic location. 

Figure 38 displays the two-dimensional ship positions, satellite ground tracks, 
and accesses for Simulation Three. In this simulation, there are four ships spaced evenly 
around the globe and transiting slowly to the East over the course of the simulation. The 
ships are labeled DDLl through DDL4, and their accesses with the satellite are shown as 
dark black traces over the satellite’s ground track. In this simulation, the ships were 
located so that most of the accesses occur over the open ocean, but there is still some 
overlapping of foreign airspace over Brazil and Indonesia. 


91 











Figure 38. Two-Dimensional Ship Positions, Satellite Ground Tracks, and Accesses 

for Simulation Three 


Figure 39 displays the Ship to Satellite accesses for Simulation Three. As 
demonstrated by the comparative size of this table to those of Simulations One and Two, 
four ship-borne lasers have the ability to generate a number of accesses to a mid¬ 
inclination satellite over a week of operation. However, it is important to note that the 
first access does not occur until July 2, thus the accesses in this scenario can be irregular 
and sometimes infrequent. With a mean duration of just over 4 minutes, these accesses 
would generally deliver power equivalent to the five minute accesses of Simulation Two. 
There is considerable variation in the length of each access with a maximum access time 
of just over 6 minutes (--389 sec) and a minimum access time of about half a minute (-30 
sec). 
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Figure 39. STK Ship to Satellite Aecess Report for Simulation Three 
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J. 


DISCUSSION 


As demonstrated by the simulations above, the eomplexities of orbital motion and 
terrestrial politieal boundaries present signifieant ehallenges to suoeessfully deploying an 
operational power beaming system. While it is ehallenging to develop and deploy an 
integrated power beaming system with dedieated ground stations and speeially designed 
satellites to take advantage of the benefits of monoehromatie illumination, an opportunity 
exists to develop this teehnology and assoeiated applieations inerementally. Assuming 
the near term deployment of high energy FELs on naval surfaee eombatants, experiments 
eould be eondueted at low power over isolated oeean regions to demonstrate the 
aequisition and traeking abilities of these vessels and the thermal effeets and potential 
power benefits for an orbiting satellite. 

The results from Simulation One demonstrate the best ease aeeess between a 
satellite at 400 km and a single ground station. In this ease, the ground station is able to 
deliver useful power levels to a sun synehronous satellite during one orbital pass every 
day. This additional power eould allow operation of the payload during eelipse or 
deerease the length of required sun soaking in the preeeding or following daylight pass, 
generally inereasing the operational use of the payload over one orbital pass. The results 
from Simulation Two spread the aeeess eoverage to two sueeessive orbital passes while 
only partially deereasing the illumination time. This would allow for more eoverage of a 
general geographie area, i.e., a war zone or disaster site requiring extra eommunieations 
or imaging. The results from Simulation Three demonstrate that multiple ground stations 
generate sporadie but generally useful aeeesses for a satellite at a moderate inelination. 
As diseussed later in Chapter VI, the same aeeesses eould be evaluated for targeting 
speeifie orbits to elear orbital debris through laser ablation or vaporization of the debris 
material. For that applieation, the length of eaeh aeeess is less important than the number 
of aeeesses, and therefore, multiple ground sites spaeed around the globe eould be a very 
useful eonfiguration for elearing orbital debris aeross a wide swath of Earth orbit. 
Eikewise, the elearanee of spaee debris is in the interest of nearly every nation, and 
therefore, it is the applieation most likely to generate the international eonsensus required 

for the development of a global power beaming eapability. 
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VI. ORBITAL DEBRIS REMOVAL 


While Chapter V investigated the use of a ground-based PEL to extend the life of 
a satellite through power beaming, Chapter VI will diseuss and evaluate the potential 
applieation of a high-powered FEE to aeeelerate the reentry of orbital debris or deerease 
the risk that they pose to operating spaeeeraft. There are four primary methods by whieh 
a laser eould affeet orbital debris. First, a laser eould be utilized to aid in the deteetion of 
non-metallie debris, whieh is diffieult to traek with radar, through illumination and 
optieal traeking. Seeond, a high-peak power laser eould ablate a small portion of the 
surfaee material, ereating a veetored veloeity ehange to lower the perigee of the orbit. 
Third, a laser eould break up the material into less massive pieees with more surfaee area. 
This method, however, generates a larger debris eloud and might only be used in lower 
orbits to ensure the break-up of objeets during reentry or to alter the reentry trajeetory to 
an uninhabited area, if possible. Fourth, the laser eould be used to heat the debris 
suffieiently to melt and then boil some of the material. As the debris material boils away, 
it should be ejeeted almost isotropieally away from the primary body, ereating a larger 
eloud of smaller debris partieles, the size of moleeules, whieh pose no risk to spaeeeraft 
and de-orbit more rapidly. All of these methods would result in an inerease in the 
atmospherie drag experieneed by the debris, and, therefore, aeeelerate the orbital deeay. 

A. BACKGROUND 

There are two types of debris in Earth orbit: meteoroids and man-made space 
debris. Meteoroids are a naturally occurring part of the interplanetary environment that 
can transit near the Earth at high velocities. Observational estimates indicate that these 
meteoroids are primarily made up of particles about 0.01 cm in diameter and that they 
move through Earth’s orbital space at an average speed of 20 km/sec [47, p. 3]. Natural 
debris is usually a “one-time” problem, i.e., the steep trajectories either enter Earth’s 
atmosphere or pass through Earth orbit and return to interplanetary space. In stark 
contrast to natural meteoroids, virtually all man-made debris, referred to here as orbital 
debris, remains trapped in Earth orbit until its final reentry into Earth’s atmosphere. 
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Orbital debris is defined by NASA as any man-made objeet in Earth orbit that no 
longer serves a useful purpose, including “derelict spacecraft and upper stages of launch 
vehicles, carriers for multiple payloads, debris intentionally released during spacecraft 
separation from its launch vehicle or during mission operations, debris created as a result 
of spacecraft or upper stage explosions or collisions, solid rocket motor effluents, and 
tiny flecks of paint released by thermal stress or small particle impacts [48].” Most of 
this orbital debris resides in high inclination orbits with velocities of 7 to 8 km/sec, and it 
can generate collisions between objects with average relative velocities of 10 km/sec 
[48]. Consequently, collisions with even small pieces of orbital debris involve 
considerable energy transfers. For comparison, in 1995, the estimated total mass of all 
the meteoroids within 2000 km of the Earth’s surface was 200 kg, and the estimated mass 
of orbital debris within the same altitudes was 2,000,000 kg [47, p. 3]. 

All orbital debris represents risk to current and future space missions, but there 
are three general categories of space debris, each with its own risk mitigation method. 
Approximately 19,000 objects larger than 10 cm are known to exist in Earth orbit. About 
9,000 of these objects are tracked in EEO and GEO orbits by the United States Space 
Surveillance Network using radar and optical observations. The orbits of these objects 
are generally well known, and some spacecraft can maneuver to reduce the chance of 
collision with them. Particles between 1 and 10 cm in diameter have an estimated 
population of about 500,000, and most of them are not tracked due to the difficulty of 
observing them with ground based radars or telescopes. This population estimate is 
generated by the intermittent radar detection of objects as small as 3mm, but ground 
based radars cannot reliably track objects of this size [48]. Without reliable track 
information, operational spacecraft cannot maneuver to avoid these mid-sized particles, 
and they are too large to be effectively shielded against. Over 10 million particles 
smaller than 1 cm in diameter are estimated to exist in Earth orbit, based on observational 
data from the external surfaces of recovered spacecraft. These objects cannot be detected 
or tracked by modem means, so spacecraft typically require shielding on some external 
surfaces to mitigate the risk of damage from collisions with them [48]. 
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Figure 40, Figure 41, and Figure 42 display an instantaneous “snap-shot” of all of 
the tracked objects in near-Earth orbits. While a small number of these objects are 
operational satellites, over 95 % of them are orbital debris. In all three images, the white 
“dots” represent the approximate instantaneous location of each tracked item, but they are 
not scaled to represent their actual size in comparison to each other or the Earth. The 
particular range and perspective of each picture was selected to adequately capture the 
orbit in question. 



Eigure 40. Objects in EEO. Erom [48] 
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Figure 41. Objects in GEO and LEO. Erom [48] 



Eigure 42. Polar View of GEO and LEO Objects. Erom [48] 
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B. 


PROBLEM SCOPE 


There are four primary types of orbital debris—^Non-Funetional Spacecraft, 
Rocket Bodies, Mission Related Debris, and Fragmentation Debris. Only one-fifth of the 
spacecraft in orbit are functional spacecraft [49, p. 21], When a GEO spacecraft reaches 
its EOL, it is usually moved to a higher or lower disposal orbit, but EEO satellites have 
traditionally been left in their orbits until they naturally decay. In both cases, these 
satellites are classified as non-functional spacecraft. Rocket bodies represent the one or 
more stages of a launch vehicle that are jettisoned as they deliver a satellite to its 
functional orbit. Usually only one rocket body is left in orbit for EEO missions, but GEO 
missions may release up to three separate rocket bodies in different staging orbits along 
the way. The presence of rocket bodies in orbit is of particular importance due to their 
typically large dimensions and the explosive potential of their residual propellants. These 
bodies generate much of the fragmentation debris in EEO. Eigure 43 is a picture of a 
Delta 2 second-stage, main-propellant tank that landed in Georgetown, Texas, in 1997. 
This tank has a mass of over 250 kg and survived reentry relatively intact [49]. 



Eigure 43. Delta 2 Second-Stage Main Propellant Tank After Reentry. Erom [48] 
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Mission-related debris is usually released in the course of a satellite’s deployment, 
activation, and operation, such as explosive bolts, spring release mechanisms, spin-up 
devices, or solid rocket slag. The amount of debris released by a typical spacecraft can 
be quite large, with one study observing the generation of 76 separate objects from a 
single Russian launch mission [49, p. 24], In Figure 44, a pellet of aluminum oxide 
(AI 2 O 3 ) solid rocket motor slag, some of the most common mission-related debris, is 
displayed. Solid rocket motors are commonly used for orbital transfer operations, but the 
relative abundance of their effluent slag has only recently been discovered and studied. 
During the burn process, large numbers of aluminum oxide particles are formed and 
ejected, but the number of particles ejected and their respective sizes are not well known. 
NASA’s Project ORION report stated that these particles are generally assumed to be less 
than 10 pm in diameter, but over one thousand can be generated in a single firing [48]. 
The particle shown in Figure 44 represents a larger and presumably less common particle 
generated during ground testing, shown for illustrative purposes. 



Figure 44. Aluminum Oxide Solid Rocket Motor Slag. From [48] 

Fragmentation debris consists of space objects created during collisions between 
other objects, explosive events, or space weapons testing. Due to the variation in ejection 
trajectories and velocities, a cloud of debris emanating from a fragmentation event 

evolves over time to spread throughout nearby orbits, resulting in a broad band of 
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material. Figure 45 displays the three phases of evolution for a eloud of orbital debris 
after a fragmentation event. In phase one, an ellipsoidal eloud of debris is formed that 
follows the general orbit of the parent body. In phase two, the veloeity and trajeetory 
variation of the debris aeeumulates to spread the ellipsoidal eloud into a torus near the 
original orbit. In phase three, the debris eloud continues to evolve into a wide band. 



Figure 45. Debris Cloud Evolution. From [50] 


Another serious source of fragmentation debris is the flaking of spacecraft paint 
under solar radiation and particle impacts. This problem was first seriously considered in 
1983, when STS-7 Space Shuttle returned to Earth with a millimeter sized crater in an 
orbiter window. Eigure 46 is a picture of this impact damage. This picture is an example 
of the risks posed by orbital debris to both manned and unmade space missions. 
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Figure 46. STS-7 Orbital Window Impact Crater due to a Paint Fleck. 

From [48] 

Ever since the launch of Sputnik in October 1957, the near-Earth environment has 
continued to accumulate orbital debris. The largest contributor to the daily growth of 
debris is the frequent fragmentation of larger objects through collision or explosive 
events. Eigure 47 displays the monthly breakdown of objects in Earth orbit by object 
type. With an international average near 75 launches each year, the number of 
spacecraft, rocket bodies, and mission related debris continues to grow, but the relative 
abundance of fragmentation debris continues to rise dramatically as collisions and 
explosions occur on-orbit. The largest fragmentation events in recent history were due to 
the Chinese Eengyun-lC (EY-IC) anti-satellite test in January 2007 and the collision 
between Iridium 33 and the derelict Cosmos 2251 in Eebruary 2009. These events 
created over 5000 fragmentation objects >10 cm in diameter and nearly doubled the 
population of cataloged orbital debris [51, p. 7]. 
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Considering the alarming rate of orbital debris generation, the era of mankind’s 
open and relatively simple aeeess to spaee may be coming to an end. Any increase of 
fragmentation events, such as through a future war with anti-satellite engagements or 
simply from the continued collisions in crowded orbits, has the potential to render those 
orbits virtually useless for generations to come. If the Chinese ASAT engagement above 
generated -3,000 pieces of debris, an anti-satellite war that destroys only 10 satellites 
could immediately double the current debris population, and this large debris field would 
spread over time to other orbits ’’near” the parent satellite. 

Currently, there are no programs for the removal of space debris from orbit, and 
the National Aeronautics and Space Administration (NASA) has only recently enacted 
guidelines to limit the creation of orbital debris. Likewise, the space debris problem will 
not “solve itself’ in the near future. The anticipated orbital lifetime of debris in the 800- 
1100 km range is on the order of 10,000 years [52, p. 576]. The space tug concept 
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discussed in Chapter V may be one method to reduee the amount of large debris, such as 
rocket bodies and non-funetional spaeeeraft, by hauling these items into lower disposal 
orbits that experienee higher atmospherie drag. Similarly, by redueing the larger parent 
objeets, mueh of the future fragmentation debris growth eould be avoided. For smaller 
debris, the most-promising, near-term method of debris removal is through the 
illumination of debris elouds with a high energy laser to lower the perigee of their orbits 
as proposed by Projeet Orion. 

C. TARGET PROPERTIES 

During the 1990s, NASA investigated the design of an orbital debris removal 
system using a ground-based laser. This study was referred to as Projeet Orion, and it 
primarily eonsidered the laser useful for the removal of orbital debris with diameters 
between 1 and 10 em. Objeets less than 1 em are extremely numerous, but they are 
diffieult to deteet. Therefore, shielding against them is eommon praetiee for all satellites. 
Objeets larger than 10 em are routinely traeked, but their numbers remain small enough 
that operational spaeeeraft ean maneuver to avoid them [53]. 

Projeet Orion foeused on five major types of mid-size debris that are most 
numerous in the orbits of eoneern. These debris types are: Na/K spheroids released from 
a leaking Russian nuelear power plant, earbon phenolie fragments, multilayered 
insulation (MLI), erumpled aluminum, and steel tank rib supports. Table 6 displays the 
typieal eharaeteristies of these debris partieles and their standard orbits. The left-most 
eolurnn lists a deseription of important orbital and physieal oharaeteristies for eaeh type 
of debris. The inelination, apogee, perigee of eaeh debris type are generalizations based 
on radar measurements of the orbit or the orbital parameters of the parent objeet. The 
approximate sizes and abundanee of eaeh debris type are based on radar or optieal 
measurements or, in some eases, theoretieal modeling. The bond albedo and optimum Cn 
are numbers that refleet the eoupling and interaetion between the laser beam and the 
objeet. The Av is a ealeulated number based on the mass and orbit of the objeets. Due to 
the age of this report, the estimated number of targets is no longer aeeurate, but the 
general eategories and orbits remain the same. As a general approximation of these 
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debris categories and the aluminum-oxide, solid-rocket-motor slag that has recently 
alarmed scientists, the calculations in this chapter will assume that the target is an 
aluminum sphere with a diameter of 5 cm. 


Debris Target Matrix 


Target 

A 

B 

C 

D 

E 



CarbM 

MU 



Deter Iplloii 

Na/K Spbere 

Pbenolle 

(Platlic/AI 

Crumpled 

steel Tank nib 



Fraonient 

Sutlacet) 

Alumimim 

Suppoil 

Inclination (dag) 

65 

87 

99 

30 

82 

Apogee (kn) 

930 

1190 

1020 

800 

1500 

Perigee Ihn) 

870 

6T0 

725 

520 

820 

A/m (cmT/gin) 

1.75 

0 7 

25 

0.37 

0.15 

Actual tiie (cm) 

1.0 

1x5 

005x30 

1x5 

1x10 

Bond albedo 

0.4 

0.02 

0.05rtl.7 

0 0510 7 

0.5 

Optiinafn Catdyne-t/JI 

6±2 

7.5*2 

5.5±2 

4i1.5 

411.5 

Av reguired (m/t) 

190 

110 

140 

90 

160 

Ettlmaled number of targala 

SOk 

20 k 

60k 

lOk 

10k 


Table 6. Debris Target Matrix. From [53] 


D, ATMOSPHERIC DRAG 

Even without human intervention, solar activity indirectly influences the orbital 
debris environment by heating the Earth’s upper atmosphere, which affects the upper 
atmospheric density and, therefore, the atmospheric drag experienced by objects in EEO. 
During the eleven-year solar cycle, fluctuations in solar intensity lead to a slow pulsation 
of Earth’s atmosphere that can accelerate the orbital decay of objects in EEO when the 
atmospheric density increases at higher altitudes. Eikewise, atmospheric drag has a 
greater effect when the orbit is at perigee, so the eccentricity of orbits tends to circularize 
over time at the same altitude as its perigee, resulting in an overall lower average altitude 
for the orbit. At lower altitudes, the object experiences greater atmospheric drag, leading 
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to further orbital decay. Therefore, any method which increases the atmospheric drag on 
an orbital debris cloud would help to decrease its orbital lifetime. 

NASA’s Project Orion sought to utilize a high-peak-power, ground-based laser to 
ablate small portions of the debris material in order to provide a change in velocity, or 
Av, that would lower the altitude of the orbit. As their moniker suggests, high-peak- 
power lasers have a very high energy per pulse, which is sufficient to cause ablation, but 
their pulses are less frequent than an FEL’s pulses. This single, high-energy pulse is 
what allows for the ablation of material and subsequent Av. 

If an FEE were utilized for orbital debris removal, it would not have enough peak- 
power to ablate the surface material of the debris. In contrast to other types of laser, an 
FEE achieves its high average power by generating frequent pulses with comparatively 
small amounts of energy per pulse. While NASA’s proposed laser for Project Orion 
generates 150 J per pulse, a 1 MW FEE will only produce about 10 mJ per pulse, but at 
~10^ times the pulse repetition frequency [53]. 

An FEE may be able to induce a faster than normal orbital decay by removing 
some of the debris material through vaporization. If enough power can be delivered in 
one orbital pass to heat some of the object enough to undergo two phase changes, from 
solid-to-liquid and then from liquid-to-gas, the remaining particle will have a smaller 
diameter and less mass than the original. As shown by Equation (IV.8), the deceleration 
due to atmospheric drag is proportional to the object’s surface area and inversely 
proportional to its mass. A less massive particle with a lower surface area may not 
experience more atmospheric drag. However, particles greater than 1 cm in size are still 
hazardous to spacecraft. Therefore, the calculations included here will determine the 
amount of energy required to melt and vaporize enough material such that the remaining 
particle is less than 1 cm. 
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E, 


CALCULATIONS 


1. Assumptions 

As in Chapters V and VII, this chapter will make a general assumption that 
atmospheric transmission for the PEL operating at 1 pm is '-'50%, so only half of the laser 
power is transmitted outside the atmosphere. These calculations will ignore turbulence 
and thermal blooming effects. Due to the small target size and pointing accuracies 
required, these effects require further study. This thesis will assume that accurate 
tracking and pointing mechanisms exist to maintain the laser beam on the targets in 
question. Improved tracking due to brighter illumination may be a significant advantage 
of the MW-class EEL. 

The calculations of this chapter approximate the target orbital debris as a 
homogeneous sphere of aluminum with a radius of 2.5 cm. This is a considerable 
oversimplification of the actual composition of most orbital debris, but aluminum is a 
fairly common material in the debris population. It is frequently utilized during 
spacecraft construction, and the aluminum oxide slag produced by solid rocket motors is 
a growing contributor to the total debris population. In terms of examining the difficulty 
of debris removal, the density, thermal conductivity, and relatively high melting and 
boiling temperatures of a metal present a challenging target for laser removal. A debris 
particle consisting of insulating material may be deflected by evaporation because they 
do not conduct heat well and the material would leave the object’s surface 
asymmetrically. 

2, STK Simulations 

In order to determine whether a high power EEE could conceivably vaporize 
orbital debris, this chapter begins with a Satellite Tool Kit (STK) simulation to determine 
the approximate access times between a hypothetical, high-altitude, ground-based laser 
and two of the prominent debris categories shown in Table 6. In this simulation, the 
altitude of the ground site does not affect the power transmitted through the atmosphere, 
but it does affect the line-of-sight to the satellite. Due to their aluminum content, relative 
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abundance, and differing orbital parameters, the materials modeled in STK were the MLI 
and the crumpled aluminum. Figure 48 displays a two-dimensional STK plot of the laser 
ground station and satellite ground tracks. The laser ground station is visible on the left, 
in the Hawaiian Islands. The high inclination orbit of the MLI debris is shown, and the 
MLI appears to the upper left of the figure. The mid-inelination orbit of the erumpled 
aluminum debris is shown near the equator, and the debris appears to the eenter-right of 
the figure. The accesses between laser and satellite tracks are shown in dark black over 
the satellites’ ground tracks. For this simulation, the lasers’ illuminations were restricted 
to times when the satellite was in sunlight and at least 10 degrees above the horizon. 



Figure 48. STK Plot of MLI and Crumpled Aluminum Debris Orbits with Laser 

Ground Station and Accesses 


Figure 49 is a summary of an STK access report for the laser ground station and 
the debris particles. In both cases, the mean duration of all accesses was over 700 sec. 
Therefore, the phase change ealeulations presented below can only be suceessful if the 
total time required is less than 700 see. 
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radiity-Hawaii*-To-Satellite-CruiDpledAl: Access Sunniary Report 


Global Statistics 


Min Duration 
Max Duration 
Mean Duration 


32 5 JUl 20>07 13144:21.893 5 Jul 2007 13:45:50.135 

6 2 Jul 2^67 01:32:^5.846 2 Jul 200*7 01 :48:05.9^4 


88.242 

930.124 

705.971 


Total Duration 


38123.413 


F.clllty-Hawali-To-S.tellite-MLC: Xccess Suamory Report 


Olobal Statistics 


Hln Duration 
Kax Duration 
Mean Duration 


19 4 Jul 2007 13:31:16.736 4 JUl 2007 13:32:46.909 
11 3 JUl 2007 03:51:34.143 3 Jul 2007 04:09:02.743 


1048.601 

709.407 

27666.862 


90.174 


Total Duration 


Figure 49. Summary of an STK Access Report for MLI and Crumpled Aluminum 

Debris With One Ground Station 


3. Calculations 

The calculations that follow attempt to evaluate the thermal properties of a debris 
particle as it is heated in daylight by solar and laser illumination. The initial temperature 
of the debris was estimated using a Microsoft Excel spreadsheet, and these calculations 
are displayed in Appendix C. All of the other calculations were performed in the 
MATLAB program, and the program used is also included in Appendix C. 

In the MATLAB code, a heat balance equation is established with inputs from 
solar and laser illumination and heat losses due to the radiation of thermal energy. The 
heat inputs due to the solar albedo and Earth’s IR radiation are not evaluated. As a first 
approximation, this model did not account for the change in propagation range as the 
debris particle passes from horizon-to-horizon or for the change in radius of the debris 
particle as material is ejected. The laser intensity is calculated at a range of 700 km and 
held constant throughout the illumination time. If the concept does not succeed without 
accounting for the losses discussed above, then it surely will not succeed when they are 
completely modeled. 
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Figure 50 displays the thermal profile of a debris particle when illuminated by a 1 
MW FEL with a beam radius at the target of 1.16 m. This laser uses similar parameters 
to the ship-home FEL discussed in Chapter V, and the beam spot size at the target is 
entirely caused by the natural diffraction of the beam. As shown below, the intensities 
produced by this laser cannot even heat the debris particle to its melting temperature 
during one orbital pass. 



Eigure 50. Debris Temperature Profile with a 1 MW Laser and a Beam Radius of 

1.16 m at the Target 

Eigure 51 displays the thermal profile of a debris particle that is illuminated by a 1 
MW EEL with a beam radius at the target of 0.25 m. This spot size could hypothetically 
be achieved using a spherical mirror that is at least 0.9 m in radius, i.e., this mirror could 
be a spherical version of the ship-borne mirror discussed above. By concentrating the 
beam to such a small area within the atmosphere, thermal blooming would have an even 
more critical effect on the beam propagation, but this simulation was intended to bound 
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the usefulness of a 1 MW laser. Even with this added intensity, the 1 MW laser still only 
heats the debris particle to about 2000 K before the heat input is matched by the radiated 
heat. 



Figure 51. Debris Temperature Profile with a 1 MW Laser and a Beam Radius of 

0.25 m at the Target 

Figure 52 displays the temperature profile of a debris particle being heated by a 
10 MW FEE with a beam radius at the target of 0.4 m. Flnlike the cases presented above, 
this laser does successfully heat the debris particle to its melting temperature and 
vaporize some, but not all, of the material that we hoped to remove. However, assuming 
the actual construction of a laser of this size, the successful atmospheric propagation of a 
laser beam with this cross-sectional area and power are questionable at best. 
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Figure 52. Debris Temperature Profile with a 10 MW Laser and a Beam Radius of 

0.4 m at the Target 


4. Discussion 

The simulations conducted above utilized many favorable approximations and 
ignored some potential loss mechanisms in the illumination of a debris particle by an 
FEL. Even without these additional losses in laser intensity at the target, the figures 
above suggest that heating large quantities of orbiting material to induce vaporization is 
not feasible for a ground-based EEE. The removal of orbital debris using a laser is 
probably best accomplished by the ablation method described by Project Orion. 
However, a ground-based PEE might be able to assist in detecting small, non-metallic 
orbital debris through illumination and optical tracking. 


112 









VII. ILLUMINATING EXTRATERRESTRIAL BODIES 


Chapters V and VI investigated the near-term engineering applications of power 
beaming and orbital debris removal for a ground-based, next-generation PEL. Chapter 
VII will discuss the application of a similar FEE to illuminate the Eunar or Martian 
surfaces for scientific study. It will also discuss the development of a larger, futuristic 
FEE for interstellar communication. Due to the high spectral density of an FEE and the 
collimation of the laser beam, it is found to be possible to match or exceed the Sun’s 
irradiance within a limited bandwidth over ranges extending out to the Martian orbit. 
Similarly, the distances involved in interstellar communications require a strongly 
collimated beam with good coherence and propagation properties. With the ability to run 
indefinitely at high-average powers, an FEE is an optimal laser to generate a repeating 
and distinctly intelligent signal for transmission to the growing number of detectable, 
“near-by” extra-solar planets that may have the ability to support life. 

A, ASSUMPTIONS 

As in Chapters V and VI, this chapter will make a general assumption that 
atmospheric transmission for the FEE operating at 1 pm is 50%, so all of the laser powers 
are halved outside the atmosphere. These calculations will ignore turbulence and thermal 
blooming effects, but due to the length of propagation outside of the atmosphere, their 
effects may seriously affect the outcome of these applications. Eikewise, some of the 
applications within this chapter assume the utilization of considerably larger lasers with 
powers that may not be feasible in the near future. Due to the futuristic nature of some of 
these applications, it may be possible that the laser could be constructed on the lunar 
surface to avoid these atmospheric losses or that some other unforeseen method of 
avoiding most atmospheric losses has been devised. A larger mirror that is already being 
considered for space exploration will use adaptive optics that can significantly reduce 
both thermal blooming and turbulence effects. 

As in previous sections, this thesis will assume that accurate tracking and pointing 


mechanisms exist to maintain the laser beam on the targets in question. As a reference 
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for the interstellar applieations in Seetion D, we can calculate the angular diameter of a 
planet-moon orbit similar to our own at interstellar ranges. The Moon orbits the Earth at 

o 

approximately 3.8 x 10 m, and the range of the probable minimum distance to intelligent 
life is about 10 light years or 10^^ m. At this range, the angular radius of the Moon’s 
orbit would be about 0.00078 arc-seconds, approximately an order of magnitude smaller 
than the 0.0072 arc-seconds angular resolution capability of the Hubble Space Telescope 
[41]. For comparison, the W. M. Keck Observatory’s Keck I telescope routinely 
achieves a pointing accuracy of 10 arc-seconds in azimuth [54]. 

For calculations regarding planetary distances and extrasolar distances, the 
distance between Earth and the other planets was based on a comparison of the semi¬ 
major axis of their orbits. This is a simplified model as the eccentricities and inclinations 
of the orbits lead to large variations in the interplanetary distances, however, it is 
assumed here that any illumination would be planned for optimal times when the orbital 
distance is small to aid the lasers propagation characteristics. 

B, SOLAR INTENSITY 

Using the Planck’s blackbody approximation given in Equation (V.2), the solar 
output within a limited bandwidth can be calculated. An FEE has a bandwidth of about 
0.1 % of its wavelength, or 1 nm for an FEE operating at 1 pm, and the Sun emits only 
about 0.07 % of its total output in this band. Starting at the surface of the Sun, which is 
assumed here to be roughly 6.96 x lO^m from the Sun’s center, this power output can be 
translated into intensity for any desired radius. While the solar energy spreads over a 
uniform sphere, a laser beam would initially remain much more collimated, until 
diffraction effects began to grow. Therefore, it may be possible for a high power laser to 
match or exceed the solar intensity within its bandwidth over some ranges. Figure 53 
displays the solar output as a function of wavelength when the sun is approximated as a 
blackbody with a temperature of 5800K. In Figure 53, the laser bandwidth is represented 
in red to show the relatively small portion of the spectrum that it occupies. 
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C. CALCULATION SETUP 


The calculations of laser beam radius (w), area (As), and intensity (!) at the target 
are based on Equations (III.l) through (III.6) and the assumptions stated above. These 
equations are restated here for easy reference. The laser beam radius (w) is given in 
terms of the distance (z), beam waist (wo), and wavelength (^) by 




(VILl) 


and the laser beam’s intensity (I) is 




P 

^ , where 


2 

= TTW 
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D, ILLUMINATING BODIES WITHIN THE SOLAR SYSTEM 


Section D will investigate the use of two notional FELs for the illumination of 
extra-terrestrial bodies within the Solar System. These FELs include the next-generation, 
ship-borne, MW-class laser that was discussed earlier and a large and futuristic laser with 
parameters beyond currently achievable technology. Due to the relatively high spectral 
intensity of these lasers, the illumination provided by them could match or exceed the 
spectral intensity of the sun inside the laser’s bandwidth. Using a specialized signal, the 
laser might even be detectable through a cross-correlation function at levels well below 
the solar, background noise. The FEE is continuously tunable so that a wide range of 
wavelengths could be used for scientific exploration or the detection of asteroids or 
meteoroids. 

Two cases will be considered here, and both of them will be assumed to operate at 
a wavelength of 1 pm with a bandwidth of 0.1 % of the wavelength. The first utilizes the 
ship-home, MW-class FEE discussed earlier that is likely to be deployed within the next 
twenty years. This laser is referred to here as Laser 1, and it produces ~1 MW of average 
power using a director mirror of about 1 m in radius. The second utilizes a similar, but 
ground-based, FEE, referred to here as Laser 2, with a larger director mirror that is 
equivalent to current optical tracking mirrors. This laser operates an average power of ~1 
MW, but its director mirror is about 10 m in radius. Although development of this laser 
would combine state of the art technology for most of the components, it could 
realistically be built in the next twenty to thirty years assuming aggressive funding. 

1. Illuminating the Moon 

o 

For these calculations, an average lunar range of 3.84 x 10 m was assumed. 
Figure 54 and Figure 55 display the laser beams’ respective radii at the lunar surface. 
Laser 1 has a radius of approximately 120 m, which could be used to illuminate specific 
topographical features of interest such as an impact crater, mountain, or canyon. Laser 2 
has a radius of approximately 16 m, which could be used for more detailed study of lunar 
topography or for power beaming to a lunar base. 
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Beam Radius vs Distance 



Figure 54. Laser 1 Beam Radius at the Moon 


Beam Radius vs Distance 



Figure 55. 


Laser 2 Beam Radius at the Moon 
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Figure 56 and Figure 57 display the laser intensities over ranges extending from 
Earth to the lunar surface. At this long range, neither laser exceeds the total AMO solar 
intensity of 1366 W/m , however, the spectral intensity of the sun is only about 96 W/m 
within the lasers’ bandwidth, which is less than an order of magnitude greater than Laser 
1 and six times less than Laser 2. 


X 10 


Intensity vs Distance 



Ligure 56. Laser 1 Intensity at the Moon 


118 












Figure 57. Laser 2 Intensity at the Moon 


In both examples, the advantages of laser illumination such as the pulse structure 
of an FEL and the ability to code the signal allows for radar-like imaging of the lunar 
surface. Likewise, in power beaming applications, the additional efficiency of solar cells 
under monochromatic illumination allows for the generation of equivalent electrical 
energy if the laser’s intensity is roughly half of the AMO solar intensity, like the 637 
W/m provided by Laser 2. Without considering atmospheric turbulence or thermal 
blooming, diffraction effects alone do not rule out the successful utilization of a laser like 
Laser 1 for scientific study of the lunar surface if a cross-correlation filter is used to 
process the return signal. Likewise, the futuristic Laser 2 produces sufficient intensity on 
the lunar surface to be an excellent candidate for beaming power to a lunar base, 
topographical mapping, or scientific illumination in the search of specific elements. 

2, Illuminating Mars 

Lollowing favorable results in lunar illumination applications, this section will 

seek to evaluate the utility of an LEL for illumination of objects as far as the Martian 

orbit. Lor this purpose, the Earth-Mars distance is approximated by the difference in the 
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semi-major axis of each orbit, about 78.3 million km, acknowledging that the Earth-Mars 
distance varies widely. It is assumed that these illuminations would be planned for times 
when the propagation distance is the smallest. In recent history, the range of closest 
approach between Earth and Mars has varied from the 55 million km approach in 2003 to 
the 101 million km approach in 1995 [55]. 

Eigure 58 and Eigure 59 display the laser beam radius for Easers 1 and 2 at Mars. 
At these ranges, the approximate radii of Easer 1 and 2 are 80 km and 8 km respectively. 
These beam sizes are large, but they are probably still useful for the illumination of large 
topographic features or specific regions, such as a potential landing zone for a spacecraft. 
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Beam Radius vs Distance 
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Figure 59. Laser 2 Beam Radius at Mars 

Figure 60 displays a comparison of the laser and total solar intensities at the 
Martian surface. As is clearly shown, even high power lasers can only compete with the 
solar spectral intensity over a relatively limited range near the Earth. At longer ranges, 
the total solar intensity exceeds Laser 1 by about 3 orders of magnitude and Laser 2 by 
about one order of magnitude. As discussed in Chapter II, the pulse structure of a MW- 
class LEL produces a peak intensity that is about 1000 times greater than the average 
intensity, therefore. Laser 1 might match, and Laser 2 could exceed, the total solar 
intensity during the macro-pulse timescale. Also, in a process common to radar and lidar 
applications, the transmitted laser pulse could contain a specific signal pattern that would 
greatly improve the detector’s signal to noise ratio (SNR) after passing through a cross¬ 
correlation filter. 
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Figure 60. Laser and Solar Intensities at Mars 

E. ILLUMINATING EXTRASOLAR BODIES 

SETI, or the Search for Extraterrestrial Intelligence, is “a human project dedicated 
to explore, understand, and explain the origin, nature, and prevalence of life in the 
universe [56].” SETI runs a number of programs searching the electromagnetic spectrum 
for signs of extraterrestrial life, and, assuming the existence of such life, there may soon 
be a day when humans become aware of advanced civilizations elsewhere in the galaxy. 
Our current understanding of science indicates that given sufficient time and the right 
conditions, life could develop on other planets. Still, mankind has not yet developed or 
deployed efficient methods of communication over such large distances. This thesis will 
avoid the political implications of communication with an unknown and not-necessarily 
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peaceful extraterrestrial civilization. This section will seek to determine what laser 
characteristics might be necessary for successful interstellar propagation of a laser 
communications signal. 

A casual survey of the night’s sky indicates that light can travel across interstellar 
distances despite considerable losses. However, the power and size of a laser built for 
communication across interstellar distances would stretch the limits of our technology. 
The development of high-average power, continuous wave FELs, when paired with a 
hypothetical, large-radius director mirror, could offer the ability to successfully propagate 
a high-powered, laser beam through the atmosphere for interstellar communications. 
While other types of laser can offer similar average powers, they tend to operate in a 
pulsed mode with a large amount of energy per pulse. This leads to large losses within 
the atmosphere, and generally, it limits their utility for transmission over interstellar 
distances. 

At extreme ranges, such as in transmission outside of the Solar System, a high- 
powered FEE can approach the spectral intensity of the Sun within its bandwidth, 
allowing for the possibility of interstellar communications. Extensions of the simulations 
presented below show that a lunar or orbiting laser of at least 10 MW with a 20 m 
director mirror will approximately match the solar spectral intensity at long ranges. 
While this is beyond current technologies, the constant development of laser, mirror, and 
space technologies indicate that a laser of this size may be achievable within the next 
century. Similarly, when combined with the previously discussed practical applications 
of high-powered lasers, it is reasonable to assume that another advanced civilization 
might already be utilizing a laser for these purposes. When beaming power for 
propulsion or communications within their own solar system, this hypothetical 
civilization may be emitting detectable signals in the direction of Earth. In this event, the 
simulations discussed here might provide some clue as to the nature of the incoming laser 
pulses. 

These simulations model the propagation of a 10 MW FEE beam over interstellar 
distances. The laser is referred to as Easer 3, and it utilizes a director mirror with a radius 
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of 10 m. This system is a “pie-in-the-sky” design to estimate the laser eharaeteristies 
required for interstellar eommunieations and to possibly reduee the thermal blooming 
effeets suffered by a smaller mirror. 

The ealeulations provide estimates of beam diameter and intensities based solely 
on diffraetion effeets. As above, only 50% of the laser power is assumed to be 
transmitted through the atmosphere, but the additional beam spreading due to thermal 
blooming is not modeled. Similarly, these models assume perfeet transmission aeross the 
interstellar medium. 

Figure 61 displays the inerease in beam radius due to diffraetion versus the range 
of transmission. The elosest star to our Sun, Alpha Centauri, is labeled at about 4 light- 
years away, and the probable minimum distanee to other intelligent life, referred to here 
as the elosest possible life, begins at least 10 light-years from Earth [57, p. 219]. For the 
10 m mirror, the beam radius inereases signifieantly over this range, growing to over 3.4 
X 10 m. This radius is almost five times larger than the radius of our Sun (-6.96 x 10 
m), but it is mueh smaller than our solar system, fitting easily inside of Mereury’s orbital 
radius (-5.8 x lO'^ m). Without knowing the eonfiguration of any potential reeeiver or 
our pointing aeouraey, it might be benefieial to produee a larger spot size to sweep out a 
larger portion of the target planetary system, but this would also further deerease the 
signal strength. Sinee habitable planets are expeeted to be within the limited orbital 
region where liquid water exists on a planet’s surfaee, this spot size may be appropriate 
for targeting a speeifie orbital band within a distant solar system. 
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Figure 61. Extrasolar Beam Radius Due To Diffraetion 


Figure 62 shows the eomparison of average laser intensity and the total solar 
intensity. The laser is about four orders of magnitude below the total intensity of our 
Sun, but as displayed in Figure 63, the eomparison of laser and solar speetral intensities 
presents a mueh more encouraging picture. At these ranges, the intensity of the Sun and 
the average laser intensity parallel each other with the laser remaining about one order of 
magnitude below the solar spectral intensity. The peak intensity of this laser should 
exceed the solar intensity by almost two orders of magnitude. With a SNR less than one, 
this laser signal will require advanced signal collection and processing to be recognized 
by a receiving civilization. The photons within these pulses should be detectable using a 
modem photon counter tube or photo-multiplier tube if gated properly to look at the right 
timing intervals. By setting a signal threshold to match the solar intensity, the receiver 
can isolate the laser transmitted photons and potentially detect the overlying signal. 
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Figure 62. Extrasolar Total Intensity Comparison 



Figure 63. Extrasolar Speetral Intensity Comparison 
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For an advanced civilization with a receiving telescope similar to the Hubble 
Space Telescope (HST), the pulsation within these limited spectral bands might be within 
their ability to detect, but this is highly dependent on the signal processing methods 
utilized. The primary mirror of the HST is 1.2 m in radius [35, p. 82]. Equation (VII. 1) 
displays a formula for the approximate number of photons received by a Hubble-like 
telescope that is 10 light years away from the Earth. Eor a wavelength of 1 pm, the 
number of photons per pulse is given by 


# photons 
pulse 


(5jd0V)(10"'" sec)(1.6jd0'‘' 




6.5jd0 


16 


(VILl) 


but these photons are spread over an extremely large beam radius. Therefore, the total 
number of photons per pulse collected by a hypothetieal receiver is 


# photons _ received 
pulse 


6.5x10*'" 


photons 

7r(l0m)^ 

pulse 

7r(3.4xl0^mj 


1 . 


(VIL2) 


A Hubble-like telescope, 10 light years away from the Earth should receive about 
one photon per pulse from the signal transmitted by the laser. While this means that the 
receiver might not detect each individual pulse, the signal could be encoded over a 
timescale of tens of microseconds, allowing the receiver to collect millions of photons at 
a time. This method could presumably be used to provide an intelligent and non-random 
pattern and communicate the evidence of our existence to other advanced civilizations. 
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F. 


DISCUSSION 


In applications ranging from power beaming to the deteetion of Earth-threatening 
asteroids, the development and deployment of a high-powered, ground-based FEE eould 
provide useful laser intensities at ranges stretehing beyond the Martian orbit. For lunar 
illumination, the beam size and intensities support many possible applieations, sueh as 
illumination for topographieal mapping, seientifie investigation of traee elements, and 
power beaming. At longer ranges near the Martian orbit, the laser intensities are less 
useful for power beaming but still offer many advantages for seientifie illumination. 
Eikewise, over the extreme ranges neeessary for interstellar eommunieations, the optieal 
quality of an FEE beam allows for good propagation and provides the potential for 
eventual eontaet with extraterrestrial life. While the lasers required for these purposes are 
well beyond our eurrent eapabilities, the potential to deteet and reeeive a similar signal 
transmitted to Earth exists and should be explored. 
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VIII. CONCLUSION 


This thesis covered a broad range of space applications for a ground or sea based 
PEL. The beam quality, high-average power, and pulse structure of an FEE give it many 
advantages over other lasers in the illumination of samples at extremely long ranges. If 
utilized for power beaming to satellites, a limited ground infrastructure of one or two 
lasers could support specific missions or orbits, such as a sun-synchronous imaging 
satellite. Eikewise, beaming power to many GEO satellites might be accomplished with 
relatively few ground stations placed in equatorial regions. However, it is concluded here 
that power beaming will not truly show its worth until satellites designed to receive laser 
illumination are constructed and deployed. By optimizing the satellite for power 
beaming, many of the complications in beam size, heating of the satellite body, and beam 
director pointing can be avoided or simplified. Due to these challenges, power beaming 
will probably first be utilized for propulsion of specialized spacecraft. The space tug 
discussed in Chapter V is a promising technology that could both reduce the costs of 
satellite launches and reduce the risks due to orbital debris by removing non-functional 
spacecraft from useful orbits. 

Orbital debris is the most challenging problem facing mankind’s future utilization 
of space. Some unmanned missions have already suffered catastrophic failures due to 
collisions with orbital debris, and the risks are mounting for future manned missions as 
the debris population grows. Recent implementation of debris reduction policies within 
NASA and throughout the international community will not reduce the current debris 
population or mitigate the growth of this population as older and larger objects continue 
to experience fragmentation events. The anticipated orbital lifetime of debris in the 800- 
1100 km range is on the order of 10,000 years, so some active removal method must be 
deployed to assure the safety of future space missions. This thesis proved that the FEE is 
not an effective option for the removal of space debris. The peak power of an FEE is too 
low to cause a Av, due to the ablation of surface material, and the average power is not 
large enough to overcome the radiative losses of thermal energy and cause vaporization 
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of the debris partiele. Still, an PEL may be useful to assist in the optieal traeking of small 
debris partieles, with sizes <10 em, that are diffieult to traek with radar alone. 

A large, ground-based FEE eould provide useful seientifie illumination of 
extraterrestrial objeets within the Solar System for ranges extending beyond the Martian 
orbit. There are many potential applieations of this teehnology ineluding the deteetion of 
asteroids, radar-like mapping of surfaee features, and the speetral seareh for speeifie 
elements. The assumed deployment of a MW-elass FEE aboard future naval eombatants 
offers a good opportunity to test the limits of this hypothetieal system and demonstrate its 
seientifie utility. Further researeh suggests that the optieal quality of an FEE beam would be 
a suitable eandidate for interstellar eommunieations. Although the direetor mirrors and lasers 
diseussed in Chapter VII slightly exeeed our eurrent teehnologies, interested Earth-bound 
researehers or organizations, sueh as SETI, might seareh for eommunieations from 
extraterrestrial sourees buried in the speetral eontent of their assoeiated stars. 

A. FUTURE WORK 

There exists a wide and interesting variety of seientifie and engineering 
ehallenges within the topies diseussed in this paper, but many of the findings presented 
here represent rough estimations of the eapabilities of a high-energy FEE. This thesis 
was intended to guide future researeh in these fields toward topies that are theoretieally 
feasible and away from those that are immediately proven impraetieal due to basie orbital 
motion or laser charaeteristies. In most oases, simplifying assumptions were made, most 
prominently in the atmospherio propagation of the laser, and future simulations or studies 
eould more aoourately measure the total effeot of the many variables listed in eaoh ease. 

In all oases, the atmosphere was assumed to remove about 50% of the laser 
power, based on the inoluded atmospherio transmission graphs in Chapter III, but 
turbulenoe and thermal blooming effeots were mostly ignored. Future researeh eould 
simulate the propagation of a high power FEE beam when transmitted vertioally through 
the atmosphere for the development of better modeling and mitigation methods for 
thermal blooming effeots. Other possible oontinuations of this researeh are diseussed 
below. 
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1 . 


Power Beaming to Satellites 


In Chapter V, this thesis diseussed various terrestrial simulations and experiments 
utilizing the nearly monoehromatic illumination of photovoltaie eells to deliver greater 
efficiencies. Future research could simulate the variation in laser intensity on the solar 
panels of an orbiting satellite. This intensity profile could then be translated into an 
actual power profile for the satellite. Similarly, as a collaboration with the NFS CubeSat 
program, future researchers could organize an experimental on-orbit illumination of a 
CubeSat by a ground based laser for measurement of the power delivered, tracking 
methods, and possibly attempt to anneal the solar cells of an older satellite. 

2, Orbital Debris Removal 

In Chapter VI, this thesis discussed the removal of orbital debris through laser 
illumination and heating. Although the simulations for this application seemed 
discouraging, significant heat can be delivered to a particle of debris. The calculations 
here restricted the debris to a pure aluminum sphere, whereas the actual debris population 
in much more complex and diverse. Depending on the specific composition of the 
targeted debris, laser illumination may be able to reduce the size of some types of orbital 
debris below the limits of current spacecraft shielding. Future research could model laser 
illumination of a more realistic debris particle to simulate the thermal and orbital effects 
or study the illumination of orbital debris to aid in detection of smaller debris particles. 

3, Illuminating Extraterrestrial Bodies 

In Chapter VII, this thesis discussed the illumination of extraterrestrial bodies 
with a high power FEL for scientific study. Future research could more accurately 
describe the desired signal and the cross correlation filter that would be required for 
radar-like applications within the solar system. The interstellar calculations were 
encouraging, but the lasers described represent a futuristic power and mirror size for the 
laser. Future research could better define the propagation effects of the interstellar 
medium and model the peak intensity of a laser for comparison to the solar intensity at 
long ranges. 
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APPENDIX A: GLOSSARY 


• Albedo - the direct reflection of solar energy from the Earth’s surface that later 
illuminates and heats a satellite in Earth orbit 

• Apogee - the point in an Earth orbit in which the orbiting object is farthest from the 
Earth 

• Argument of Perigee - the angle between the line extending from the center of the 
Earth to the ascending node of an orbit and the line extending from the center of the 
Earth to the perigee point in the orbit measured from the ascending node in the 
direction of motion of the satellite 

• Attitude - the orientation of a spacecraft as determined by the inclination of its axis 
with respect to a fixed reference point on Earth 

• Ascending Node - the point on the ground track of a satellite’s orbit where the 
subsatellite point (SSP) crosses the equator passing from the Southern to the Northern 
Hemisphere 

• Atmospheric Turbulence - refers to fluctuations in the index of refraction resulting 
from small temperature fluctuations 

• Beam Waist - the smallest spot size along the propagation path of a collimated or 
convergent beam where the on -axis irradiance is a maximum and the phase front 
radius of curvature is infinite 

• Beam Wander - deflection of the beam centroid from the boresight of the optical 
wave 

• Coherence (spatial) - the ability of a light beam to interfere with a spatially shifted 
version of itself 

• Coherence (temporal) - the ability of a light beam to interfere with a delayed version 
of itself 

• Collimated Beam - a Gaussian -beam wave with phase front radius of curvature that 
is infinite at the transmitter (beam waist is located at the transmitter) 

• Delta -V - the change in velocity vector of an object 

• Descending Node - the point on the ground track of a satellite’s orbit where the 
subsatellite point (SSP) crosses the equator passing from the Northern to the Southern 
Hemisphere 

• Downlink - the communications link from a satellite down to one or more ground 
stations 

• Eccentricity - a fixed constant for each type of conic section representing an orbit’s 
“roundness” 

• Ecliptic - the apparent path that the Sun traces out in the sky during the year 
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Figure 64. Electromagnetic Spectrum. From [58] 

• Ephemeris - a table of values that gives the positions of astronomical objects in the 
sky at a given time 

• Epoch - specifies the reference time for a set of orbital elements 

• Elliptical Orbit - a Keplerian orbit with an eccentricity of less than one; this includes 
the special case of a circular orbit, with an eccentricity of zero 

• Far Field - also called the Frau nh ofer Region. 

• Fraunhofer Region - the region outside the near -field region where the angular 
field distribution is essentially independent of distance from the source 

• Fresnel Region - the close -in region of an antenna where the angular field 
distribution is dependent upon the distance from the antenna 

• Gaussian Beam - a beam of electromagnetic radiation whose wave front is 
approximately Gaussian (parabolic) at any point along the path and whose transverse 
field intensity over any wavefront is a Gaussian function of distance from the beam 
axis (refers to the lowest order mode) 

• Geostationary Orbit - an orbit at zero degrees inclination at an altitude of 35,687 km 
that precisely matches the Earth’s angular velocity, thereby keeping the satellite at a 
stationary point in the sky when viewed from the ground 

• Ground Station - a terrestrial terminal station designed for communication with 
spacecraft 
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• Ground Track - the path on the surfaee of the Earth traeed out direetly below a 
satellite 

• Heterodyning - mixing or multiplying two oscillating waveforms, used in 
modulating or demodulating signals 

• Inclination - the angle the orbit plane makes with the equatorial plane 

• Index of Refraction - the ratio of the phase velocity of light in a vacuum to that in a 
specified medium 

• Intensity - a measure of the power per unit area of an illumination source 

• Isp - specific impulse; a way to describe the efficiency of jet or rocket engines. It 
represents the change in momentum per unit amount of propellant used 

• Keplerian Elements - seven numbers used to describe the orbit of an object. These 
numbers define an ellipse, orient it about the Earth, and place the satellite on the 
ellipse at a particular time. The basic orbital elements are: epoch, orbital inclination, 
RAAN, argument of perigee, eccentricity, mean motion, and mean anomaly 

• Lorentz Force - the force on a point charge due to electromagnetic fields, given by 
the equation F = qi^E -\-vxB^ 

• Near -Field Region - also called the Fresnel Region 

• Orbital Debris - any man -made object in Earth orbit which no longer serves a useful 
purpose 

• Orbital Elements - see Keplerian elements 

• Perigee - the point in the orbit that is nearest the Earth 

• Perturbation - deviation from a normal, idealized, or undisturbed motion 

• Phase front - a surface of constant phase (also called a wave front) 

• Plane Wave - a wave in which the phase fronts form parallel planes 

• Rayleigh Range (zr) - defines the distance along the propagation axis on either side 
of the beam waist up to the point where the beam spot size doubles 

• Refraction - a change in the direction of propagation of any wave 

• Resolution - the ability of an optical system to resolve details in an object being 
imaged 

• Right Ascension of the Ascending Node (RAAN) - the angle between the line 
extending from the center of the Earth to the ascending node of an orbit and the line 
extending from the center of the Earth to the vernal equinox measured from the vernal 
equinox eastward in the Earth’s equatorial plane 

• Scintillation - refers to either temporal or spatial fluctuations in the irradiance of an 
optical wave (e.g., star twinkle) 
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• Spontaneous Emission - the process by which an electron, initially in a higher 
energy state, spontaneously drops to a lower energy state by emitting a photon whose 
energy exactly matches the difference in states 

• Spot Size - refers to the 1/e beam radius of a Gaussian beam 

• Stimulated Emission - the process by which an electron, stimulated in the presence 
of an electromagnetic wave, emits a photon with the same phase, frequency, 
polarization, and direction of travel as the original wave 

• Sub -Satellite Point - the point on the Earth’s surface directly below a satellite, i.e., 
on a line between the center of the Earth and the satellite center 

• Sun Soak/Bathe - When a satellite orients its solar panels towards the sun for 
maximum battery charging. Often, mission requirements might require prolonged 
periods of operation without the solar arrays optimally aligned, normal to the incident 
solar rays. 

• Sun -Synchronous Orbit - a satellite orbit that matches the precession of perigee to 
counter the daily orbital motion of the Earth around the Sun; this produces an orbit 
that passes overhead at the same local time 

• Synchrotron Radiation - electromagnetic radiation generated by the acceleration of 
relativistic charged particles by a magnetic field 

• Vernal Equinox - the direction of the Sun in space when it passes from the southern 
hemisphere to the northern hemisphere and appears to cross the Earth’s equator 
(about 20 March). The vernal equinox is the reference point for measuring angular 
distance along the Earth’s equatorial plane (right ascension), and one of two angles 
used to locate objects in orbit the other is declination) 

• Wave front - see phase front 
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APPENDIX B: SIMULATION SETTINGS & RESULTS 


This appendix outlines the specific settings within STK version 8 that were 
utilized in the simulations of Chapter V. All simulations were conducted with a 
Coordinate Epoch of 1 Jan 2000 11:58:55.816 UTCG for each scenario, and the scenarios 
ran over a one week period from 1 Jul 2007 to 8 Jul 2007. 



Figure 65. Sun Synchronous Satellite from Simulation One 
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Figure 67. Global In-View Azimuth, Elevation, and Range Summary for Simulation 

One 


Propagator: |j2Perturbation 




Central Body: Earth 


Start rime: 
Stop Time: 
Step Size: 

Orbit Epoch: 
Coord Epoch: 
Coord Type: 
Coord System: 


IJJ 2007 12:OC:00.000 UTCG ^ 


8 JJ 2007 12:OC:00.000 UTCG m 


GO sec 


1 JJ 2007 00:OC:00.000 UTCG ™ 


Senimajor Axis 




G939.13km 


1 Jan 2G\j 0 11:58:55.816 UTCG mj |Eccentricity 

3 


Qassical 


J2000 


"3 


Argument cf Perigee |o deg 


Prop Specific: Special Options... | 


|RA\N 

|Trus Anomaly ^ |0deg 


Initial 


|2.26667e-018 ^ 

Inclination |97.6402 deg rs?| 


y I |2?8.574 deg 


Figure 68. Sun Synchronous Satellite from Simulation Two 
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Figure 69. Faser Equipped Ship in Simulation Two 
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Figure 70. Global In-View Azimuth, Elevation, and Range Summary for Simulation 

Two 
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Figure 71. Mid-Inclination Satellite with Zero RAAN from Simulation Three 
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Figure 72. Faser Equipped Ship DDFl from Simulation Three 
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Figure 73. Laser Equipped Ship DDL2 from Simulation Three 
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Figure 74. Laser Equipped Ship DDLS from Simulation Three 
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Figure 75. Laser Equipped DDL4 from Simulation Three 
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Figure 76. Global In-View Azimuth, Elevation, and Range Summary for Simulation 

Three 
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APPENDIX C: ORBITAL DEBRIS THERMAL ESTIMATE 


Figure 77 is a Microsoft Excel spreadsheet used to approximate the initial 
temperature of the aluminum sphere in Chapter VI. 
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Figure 77. Thermal Temperature Calculations for a 5 cm Sphere of Aluminum at an 

Altitude of 700 km 
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The following MATLAB code was used to calculate the temperature of a debris 
particle illuminated during the daylight. The results are discussed and included in 
Chapter VI. 

% LCDR Benjamin Wilder 

% Master's Thesis Work- Time to vaporize a debris particle 
% Summary- This MATLAB m file was created in MATLAB 7.7.0 (R2008b). It 

% assumes a homogeneous, spherical debris particle of pure aluminum. 


clear all; close all; clc; 
% Constants 


To = 

270; 

"6 

da = 

2.702; 

"6 

d = 

da*(100*3)/lOOO; 

"6 

cs = 

0.910*1000; 

"6 

Tm = 

933; 

"6 

Lf = 

397*1000; 

"6 

Tb = 

2723; 

"6 

Lv = 

11400*1000; 

"6 

sigma = 5.67e-8; 

"6 

e = 

0.7; 

"6 

% Inputs 


r = 

2.5/100; 

"6 

rf = 

1/100; 

"6 

Pltl 

= 5e6; 

"6 

rbt 1 

o 

II 

"6 

% Calculations 


Vo = 

4*pi*(r*3)/3; 

"6 

Ao = 

pi*r*2; 

"6 

As = 

4*pi*r*2; 

"6 

m = 

Vo*d; 

"6 

Vf = 

4*pi*(rf*3)/3; 

"6 

mf = 

Vf*d; 

"6 


Initial Temp in K 

Density of aluminum in g/cm''3 

Density of aluminum in kg/m''3 

Specific Heat Capacity in J/K-kg 

Melting Temp in K 

Latent Heat of Fusion in J/kg 

Boiling Temp in K 

Latent Heat of Vaporization 

Boltzman's constant in W/m^2-K^4 

emissivity of aluminum 


Initial radius of the debris in m 
Final radius of the debris in m 
Power of the laser at the target 
Beam radius at the target in m 


Initial volume of the debris in m''3 
Cross sectional area of the debris 
Surface area of the particle 
Mass of the particle in kg 
Final volume of the debris in m^3 
Final mass of the particle in kg 


Ilaserl = Pltl/ (pi*rbtl''2) ; % Laser 1 intensity at the target in W/m''2 


Pinl = Ilaserl*Ao+1366*Ao; % Power received by the debris particel in W 


% Calculate the time to heat to melting temperature 


T(l) = To; 
t(l) = 0; 
deltat=0.01; 
n = 1; 

while T(n) < Tm 

n = n + 1; 

t(n) = t(n-l) + deltat; 


Initialize temp 

Initialize time in seconds 

Time step in seconds 

Creates an array number for t and T 


% Current time in seconds 
Prad = e*sigma*((T(n-1)-3)^4)*As; % Radiated power by the particle 
Pdif = (Pinl-Prad); % Power difference 

delta! = Pdif*deltat/(m*cs); % Change in temp due to power (in K) 

T(n) = T(n-l) + delta!; % Current temp 

if t(n) > 700 
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break 

end 

end 
T(n) 
t (n) 

% Calculate the time to complete the 
Q = 0; 

Qm = m*Lf; 
while Q < Qm 

if t(n) > 700 
break 

end 

n = n +1; 

t (n) = t(n-l) + deltat; 

Prad = e*sigma*((T(n-1)-3)^4)*As 
Pdif = (Pinl-Prad); 

Q1 = Pdif*deltat; % 

Q = Q + Ql; 

T(n) = Tm; 

end 


phase transition (melting) 

% Initializing the Q counter 
% Heat required to melt in J 


% Time counter in seconds 
% Radiated power by the particle 
% Power difference 

Change in temp due to power (in K) 
% Heat counter in J 
% Keeps track of T for t 


% Calculate the time to heat to boiling temperature 
while T(n) < Tb 

n = n + 1; % Counter 

t(n) = t(n-l) + deltat; % Current time in seconds 

Prad = e*sigma*((T(n-1)-3)^4)*As; % Radiated power by the particle 
Pdif = (Pinl-Prad); % Power difference 

delta! = Pdif*deltat/(m*cs); % Change in temp due to power (in K) 

T(n) = T(n-l) + delta!; % Current temp 

if t (n) > 700 
break 

end 

end 


% Calculate the time to complete the phase transition (boiling) 

Q2 = 0; % Initializing the Q counter 

Qb = m*Lv; % Heat required to melt in J 

while Q2 < Qb 

if t(n) > 700 

break % ends the loop if program is too long 

end 

n = n + 1; 

t(n) = t(n-l) + deltat; % Time counter in seconds 

!(n) = Tb; 

Prad = e*sigma*((T(n-1)-3)^4)*As; % Radiated power by the particle 
Pdif = (Pinl-Prad); % Power difference 

Qdif = Pdif*deltat; % Change in temp due to power (in K) 

Q2 = Q2 + Qdif; % Heat counter in J 

End 

% Plot the results 
figure(l); plot(t,T); 

title (' Temperature vs Time'); ylabel('Temp (K)' ); xlabel('Time (s)'); 
axis([0 700 270 Tb+100]); Tmline(l); Tbline(l); Tbtext(l); Tmtext(l); 
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APPENDIX D: ILLUMINATION MATLAB CODE 


The following code was created in support of the calculations discussed in 
Chapter VII. The contributions of LCDR Sean Niles and LT Justin Jimenez in 
programming advice and debugging are greatly appreciated. [59] [60] 

% LCDR Benjamin Wilder, USN 

% Masters thesis work- Illumination of extraterrestrial bodies 
% This work was conducted within MATLAB 7.7.0 (R2008b) 

% Summary- This m file will utilize Planck's analytical expression for 
% the Stefan-Boltzman Law to calculate the irradiance of the Sun by 
% approximating it as a blackbody with a temperature of 5800 K. 

% Specifically, the power is calculated within the limited 1 nm 
bandwidth 

% of an FEL operating at 1 micrometer. This irradiance is plotted 
versus 

% radius from the Sun and compared to the irradiance of free elctron 
% lasers of varying mirror size and total power 

% The contributions of LCDR Sean Niles and LT Justin Jimenez in 
programming 

% advice and debugging is greatly appreciated. 

clear all; close all; clc; format long e; 

%% Constants and Parameters 
% FEL paramters 

wol =1; % Mirror radius in m 


wo2 = 10; 
wo3 = 20; 

Pi = 0.5e6; 

atmosphere 

P13 = 5e6; 

wavelength = le-6; 

zol = pi*wol^2/wavelength 

zo2 = pi*wo2''2/wavelength; 

zo3 = pi*wo3^2/wavelength; 

bandwidth = 0.1; 

% Given or Assumed Constants 
h = 6.6260693e-34; 
k = 1.380658e-23; 
c = 2.9978e8; 

T = 5800; 


% wavelength in m 
% Rayleigh Range 


% Power of the Laser outside the 


% bandwidth in % 


% Planck's Constant in J-sec 
% Boltzman's cosntant in J/K 
% Speed of light in m/sec 
% Blackbody Temperature in Kelvin 


9- 9- 
0 o 


Create the total blackbody curve 


% wavelength range 

lambda = linspace(0.5e-7,5e-6, 1000) ; 


% compute intensity per m at each point 
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M = (2*pi*c"2*h)./(lambda.^5).*1./(exp((h*c)./(lambda.*k*T))-1); 

% integrate area under curve 
I_func = 

@ (lambda) (2*pi*c^2*h) ./ (lambda .^5) .*l./(exp((h*c) ./ (lambda . *k; . . . 
*T))-l); 

I_band = quad(I_func,0.9995e-6,1.0005e-6) 

I_tot = quad(I_func,0.5e-6,lle-6) 

Percentage = I_band*100/I_tot; 


%% Calculate the Sun's total power 
laser 

% Constants from Orbital Mechanics 
Rs = 6.96e8; % 

Re = 149.6e9; % 

Rmoon = 384.4e6; % 

Rmars = 227.9e9; % 

Rjupt = 778.6e9; % 

Rsatn = 1.433el2; % 

Ruran = 2.872el2; % 

Rnept = 4.495el2; % 

Rpluto = 5.87el2; % 

Psun = 4e26; % 

output 

Isunsurf = Psun/(4*pi* ( (Rs)''2) ) % 

surface 

%% Creating the laser and solar intensity comparison 

FF = 8001; % Number of steps to be calc'd 

Minradius = 1.4ell; 

Maxradius = 3ell; 

Diffradius = (Maxradius-Minradius)/(FF-1); 
for F=1:FF 

r = Minradius + (F-1)*Diffradius; 

R(F) = r; 

Isuntot(F) = Isunsurf*(Rs/r)^2; % Total solar intensity at r in 

W/m-'2 

Iband(F) = I_band*(Rs/r)^2; % Spectral solar intensity at r in 

W/m''2 

if r>= Re 

wl = wol*sqrt(1+((r-Re)/zol)^2) ; 

Ilaser(F) = PI/(pi* (wl)''2) ; % Laser Intensity at radius r 

w2 = wo2*sqrt(1+((r-Re)/zo2)^2) ; 

Ilaser2(F) = PI/(pi*(w2)*2); % Laser Intensity at radius r 

w3 = wo3*sqrt(1+((r-Re)/zo3)^2); 

Ilaser3(F) = P13/(pi*(w3)^2); % Laser Intensity at radius r 


% The percentage of solar output 
% within the laser bandwidth 

for later comparison output to the 


Textbook 

Solar Radius in m 
Earth Orbital Semi-major Axis 
Lunar Semi-major Axis 
Mars Orbital Semi-major Axis 
Jupiter Orbital Semi-major Axis 
Saturn Orbital Semi-major Axis 
Uranus Orbital Semi-major Axis 
Neptune Orbital Semi-major Axis 
Pluto Orbital Semi-major Axis 

The estimated total solar power 
The Total Solar Intensity at the 


else 

end 

end 
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%% Plot the results 
figure(1); 

plot(lambda,M, '-b' , 'linewidth' ,2); 

xlabel( 'Wavelength (m) ' ) ; 
ylabel ( 'M (W/m''3) ' ) ; 

title ('Solar Blackbody Approximation'); 
hold on; laserbandwidth(1); 

annotation (' textbox' , [0.3 0.4 0.3 0.1],... 

'Interpreter' , 'tex' , 'Edgecolor' , 'none' , ... 

'String', ['% of solar output within the laser bandwidth is ' ... 

num2str(Percentage)], 'FitBoxToText' , 'on' ); 

figure(2); % plot total solar intensity vs 

laser 

loglog(R, Isuntot, '-k' , R, Ilaser, '.r', R, Ilaser2, '*b' , Re, leS, 

'ob' , Rmars, leS, 'ob' ); 

h = legend (' Solar' , '1 MW, 1 m' , '1 MW, 10 m' ,1); earthbox(2); 
marsbox(2); 
hold on; 

title ('Laser vs Total Solar Intensity Comparison'); 
ylabel( 'Intensity (W/m^2)' ); xlabel( 'Radius (m)' ); 

figure(3); 

loglog(R, Iband, '-k' , R, Ilaser, '.r', R, Ilaser2, '*b' , Re, leS, 'ob' , 
Rmars, leS, ' ob' ); 

h = legend( 'Solar' , '1 MW, 1 m' , '1 MW, 10 m' ,1); earthbox(3); 
marsbox(3); 
hold on; 

title ('Laser vs Solar Spectral Intensity Comparison'); 
ylabel (' Intensity (W/m''2)'); xlabel (' Radius (m) ' ) ; 


%% Compute and plot the Spot Size 
0=8001; 

Rmin=le3; 

Rmax=Rmars; 

Diffradius = (Rmax-Rmin)/(0-1) ; 
for F=1:0 

r = Rmin + (F-1)*Diffradius; 

R1(F) = r; 

wl (F) = wol*sqrt (1+(r/zol)''2 ) ; 
A1(F) = (pi* (wl (F) ) ■'2) ; 

Ilaser(F) = Pl/Al(F); 
r 

w2(F) = wo2*sqrt(1+(r/zo2)*2); 
A2(F) = (pi*(w2(F))*2); 

Ilaser2(F) = P1/A2(F); 
r 

w3 (F) = wo3*sqrt (1+(r/zo3)''2) ; 
end 


versus distance 


% Beam radius in m 
% Spot size in m''2 

% Laser Intensity at radius 

% Beam radius in m 
% Spot size in m''2 

% Laser Intensity at radius 

% Beam radius in m 


% plot the results 

figure(4); 

plot(R1,wl, '-b' ); 
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h = legend ('1 MW, 1 m',2); 
set(h, 'Interpreter' , 'none' ); 
xlabel( 'Distance (m)' ); 
ylabel( 'Beam Radius (m) ' ) ; 
title ('Beam Radius vs Distance') 

figure(5); 

plot(R1,A1, '-b' ); hold on; 
h = legend (' 1 MW, 1 m',2); 
set(h, 'Interpreter' , 'none' ); 
xlabel( 'Distance (m)' ); 
ylabel ('Spot Size (m^2)'); 
title ('Spot Size vs Distance'); 

figure ( 6); 

plot(Rl,Ilaser, '-b' ); hold on; 
h = legend ('1 MW, 1 m',2); 
set(h, 'Interpreter' , 'none' ) ; 
xlabel( 'Distance (m)' ); 
ylabel( 'Intensity (W/m^2) ' ) ; 
title (' Intensity vs Distance'); 

figure(7); 

plot(Rl,w2, '-r' ); 

h = legend (' 1 MW, 10 m',2); 

set(h, 'Interpreter' , 'none' ); 

xlabel( 'Distance (m)' ); 

ylabel( 'Beam Radius (m) ' ) ; 

title ('Beam Radius vs Distance') 

figure (8); 

plot(Rl,A2, '-r' ); hold on; 
h = legend ('1 MW, 10 m',2); 
set(h, 'Interpreter ' , 'none ' ); 
xlabel( 'Distance (m)' ); 
ylabel ('Spot Size (m^2)'); 
title('Spot Size vs Distance'); 

figure(9); 

plot(Rl,Ilaser2, '-r' ); hold on; 
h = legend (' 1 MW, 10 m',2); 
set(h, 'Interpreter' , 'none' ); 
xlabel( 'Distance (m)' ); 
ylabel ( ' Intensity (W/m''2) ' ) ; 
title (' Intensity vs Distance'); 


The following m file supports the Extrasolar illumination discussions in Chapter 
VII, Section E. 

% LCDR Benjamin Wilder, USN 

% Masters thesis work- Illumination of extraSOLAR bodies 
% This work was conducted within MATLAB 7.7.0 (R2008b) 

% Summary- This m file will compare the spectral intensities of the 
laser 

% and the Sun within the laser bandwidth at extreme ranges. Further 
plots 

% will show the laser spot size versus range from the Earth at star 
systems 

% of interest. 

clear all; close all; clc; format long e; 

%% Constants and Parameters 
% FEL paramters 

wol = 10; % Futuristic Mirror radius in m 

wo2 = 100; % Pie-in-the-sky Mirror radius in m 

Pi = 5e6; % Power of the Laser outside the atmosphere 


wavelength = le-6; 

"6 

wavelength in m 

zol = pi*wol^2/wavelength; 
zo2 = pi*wo2''2/wavelength; 

"6 

Rayleigh Range 

bandwidth = 0.1; 

% Given or Assumed Constants 

"6 

bandwidth in % 

h = 6.6260693e-34; 

"6 

Planck's Constant in J-sec 

k = 1.380658e-23; 

"6 

Boltzman's cosntant in J/K 

c = 2.9978e8; 

"6 

Speed of light in m/sec 

T = 5800; 

"6 

Blackbody Temperature in Kelvin 


%% Create the total blackbody curve 
% wavelength range 

lambda = linspace(0.5e-7,5e-6, 1000) ; 

% compute intensity per m at each point 

M = (2*pi*c"2*h)./(lambda.^5).*1./(exp((h*c)./(lambda.*k*T))-1); 

% integrate area under curve 
I_func = 

@(lambda)(2*pi*c^2*h)./(lambda.^5).*l./(exp((h*c)./(lambda.*k ... 
*T))-l); 

I_band = quad(I_func,0.9995e-6,1.0005e-6) 

I_tot = quad(I_func,0.5e-6, lle-6) 


%% Calculate the Sun's 
laser 

Rs = 6.96e8; 

Re = 149.6e9; 

Rpluto = 5.87el2; 
Rpluto = Rpluto-Re; 
Ralphacentari=4; 


total power for later comparison output to the 

% Solar Radius in m 
% Earth Orbital Semi-major Axis 
% Pluto Orbital Semi-major Axis 
% Shortest range to pluto from earth in m 
% Range to alpha centari in light years 
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Rlife= 10; 


% Probable distance to nearest intelligent 
% life in light years 
Rac = Ralphacentari*9.4605284el5; % Rac in m 
Rlife= Rlife*9.4605284el5; % Rlife in meters 


Psun = 4e26 % The estimated total solar power output 

Isunsurf = Psun/(4*pi* ( (Rs)''2) ) % Total Solar Intensity at the surface 
%% Creating the laser and solar intensity comparison 


FF = 8001; % Number of steps to be calc'd 

Diffradius = (Rlife*l.1-Re)/(FF-1); 
for F=1:FF 

r = Re + (F-1)*Diffradius; 

R(F) = r; 


Isuntot(F) = Isunsurf*(Rs/r)^2; % Total solar intensity 

Iband(F) = I_band* (Rs/r)''2; % Spectral solar intensity 

wl(F) = wol*sqrt(1+(r/zol)^2); % Beam radius in m 

Ilaser(F) = PI/(pi* (wl (F) )''2) ; % Laser Intensity at 

w2(F) = wo2*sqrt(1+(r/zo2)*2); % Beam radius in m 

Ilaser2(F) = PI/(pi*(w2(F))^2); % Laser Intensity at 


at r in W/m^2 
at r in W/m''2 

radius r 

radius r 


end 


%% Plot the results 

% Plot the laser beam radius versus distance 
figure(1); 

plot(R,wl, '-b' , Rac, 2e9, Rlife, 2e9); 

createlinel(1) ; createline2(1) ; actextl(l); cpitextl(l); 
mirr1(1); 

xlabel( 'Distance (m)' ); 

ylabel( 'Beam Radius (m) ' ) ; 

title ('Beam Radius versus Distance'); 


% plot the laser intensity versus the total solar intensity 
figure (2); 

loglog(R, Isuntot, '-k' , R, Ilaser, '.r'); 

hold on; plutotext2 (2); plutoline2 (2); actext2(2); acline2 (2); . . . 

cpitext2 (2); cpiline2 (2); solar2 (2); mirrorl (2); 
title ('Laser vs Total Solar Intensity Comparison'); 
ylabel( 'Radius (m) ' ) ; xlabel( 'Distance (m) ' ) ; 

% plot the laser intensity versus the solar intensity in the laser 
% bandwidth 
figure(3); 

loglog(R, Iband, '-k' , R, Ilaser, '.r'); 

hold on; plutolineS(3) ; acline3(3); solar3(3); 

cpiline3(3); plutotext3(3); cpitext3(3); actext3(3); mirror3(3); 
title ('Laser vs Solar Spectral Intensity Comparison'); 
ylabel (' Intensity (W/m^2)'); xlabel (' Distance (m)' ); 
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